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Since  1950’ s,  many  metallurgical  researches  on  fatigue  crack  initiation  were  undertaken.  Now,  the 
difficulty  is  to  see  these  multitudinous  facts  as  a  related  and  connected  whole,  largely  because  of  the  lack  of 
a  sound  general  theory.  This  present  research  is  an  attempt  to  develop  such  a  general  theory.  It  is  essential 
in  estimating  (interpolating  and  extrapolating)  fatigue  lives  under  different  design  conditions.  Fatigue 
failure  has  been  observed  at  temperatures  as  low  as  l.T^K.  This  indicates  that  although  corrosion  and 
diffusion  may  have  important  influence  on  but  not  necessary  to  fatigue  failure.  This  leaves  the  local  stress 
and  strain  to  provide  the  basic  mechanism  of  fatigue. 

Based  on  the  hint  of  intrusions  and  extrusions  of  fatigue  specimens,  a  micromechanic  theory  of 
high-cycle  fatigue  crack  initiation  has  been  developed.  This  theory  has  been  found  to  have  extensive 
metallurgical  verifications.  This  theory  has  been  applied  to  study  the  interaction  of  multiple  fatigue  bands 
and  the  propagation  of  fatigue  bands  across  die  grain  boundary.  This  propagation  of  fatigue  bands  has 
been  found  to  depend  greatly  on  the  misorientation  of  the  nei^boring  crystal.  The  effect  of  the  elastic 
anisotropy  of  the  component  crystals  of  a  polycrystal  is  analyzed. 
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Abstract 


Since  1950's,  many  metallurgical  researches  on  fatigue  crack  initiation  were  undertaken. 
Now,  the  difficulty  is  to  see  these  multitudinous  facts  as  a  related  and  connected  whole,  largely 
because  of  the  lack  of  a  sound  general  theory.  This  present  research  is  an  attempt  to  develop  such 
a  general  theory.  It  is  essential  in  estimating  (interpolating  and  extrapolating)  fatigue  lives  xmder 
different  design  conditions.  Fatigue  failure  has  been  observed  at  temperatures  as  low  as  l.T^K. 
this  indicates  that  although  corrosion  and  diffusion  may  have  important  influence  on  but  not 
necessary  to  fatigue  failure.  This  leaves  the  local  stress  and  strain  to  provide  the  basic  mechanism 
of  fatigue. 

Based  on  the  hint  of  intrusions  and  extrusions  of  fatigue  specimens,  a  micromechanic 
theory  of  high-cycle  fatigue  crack  initiation  has  been  developed.  This  theory  has  been  found  to 
have  extensive  metallurgical  verifications.  This  theory  has  been  applied  to  study  the  interaction  of 
multiple  fatigue  bands  and  the  propagation  of  fatigue  bands  across  the  grain  boundary.  This 
propagation  of  fatigue  bands  has  been  found  to  depend  greatly  on  the  misorientation  of  the 
neighboring  crystal.  The  effect  of  the  elastic  anisotropy  of  the  component  crystals  of  a  polycrystal 
is  analyzed. 
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I.  Introduction 


Stresses  in  polycrystal  under  thermal  or  mechanical  loadings  are  highly  heterogeneous. 
These  stresses  causes  localized  plastic  deformation,  which  often  appears  on  the  free  surface  as  slip 
bands.  These  bands  developed  under  monotonic  loadings  are  known  as  shear  band  and  those 
developed  under  cyclic  loadings  are  known  as  fatigue  bands.  This  analysis  of  slip  bands,  is 
applicable  to  those  developed  under  both  thermal  and  mechanical  loadings. 

In  Section  II,  the  mechanics  of  shear  band  formation  is  given.  In  Section  III,  the 
micromechanism  of  fatigue  band  is  shown.  This  gives  a  physical  model  of  high-cycle  fatigue 
(HCF)  crack  initiation  (Lin,  1992).  The  interaction  of  multiple  fatigue  bands  in  the  most  favorably 
oriented  crystal  at  a  free  surface  is  shown  in  Section  IV  (Lin  and  Chen,  1989).  Propagation  of 
fatigue  bands  across  grain  boundaries  is  given  in  Section  V.  HCF  crack  initiation  of  intermetallic 
compounds  is  shown  in  Section  VI.  Anisotropy  of  elastic  constants  of  crystal  of  aluminum  and  its 
alloys  is  small  and  was  neglected.  But  this  anisotropy  of  other  metals  like  titanium  and  NisAl 
intermetallic  alloys  are  not  small.  Hence,  the  effect  of  this  anisotropy  on  jwlycrystal  fatigue  bands 
is  analyzed  in  Section  VII. 


11.  Shear  Band  Formation 

A  long  thin  slice  of  metal  of  uniform  rectangular  cross-section  is  shown  in  Figure  1, 
embedded  in  an  infinite  isotropic  elastic  medium.  The  shear  stress  field  xi2  caused  by  a  uniform 
plastic  shear  strain  ejj  has  been  calculated  (Lin,  1972, 1992). 


1^2 


Rgure  1 

Thin  slice  cross-section. 
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and  is  given  (with  the  subscript  "12"  dropped)  as 


2pe^>v 

(Xj  +  d) 

(X2 -w^)-(Xi+d)^ 

1 

3l(l-  V) 

(j^  +  d)^  +  (Xj 

-  wf’ 

(Xj  +  d)^  +  (Xj 

,  +  w)2] 

ixi  -  d)  {Xi -df  -  (Ji2  - 
^  [(jTi  -  df  +  (JC2  +  -  df  +  {xl  -  w)^ 


(1) 


Let  e^w  be  kept  constant,  let  increase  and  let  w  decrease  to  zero,  then 


2pc^w 

1 

1 

^  1 
1 

1  (Xi  +  d)[(Xi  +  d)^-X2] 

Jt(l-  v) 

[{Xi-df  +xlj 

[(Xj  +  d)^  +  X2] 

The  resolved  shear  stress  due  to  the  presence  of  an  edge  dislocation  along  the  ;c3-axis  (0,0,X3)  with 
a  Burger's  vector  pointing  in  the  xi-direction,  is  given  by  Hirth  and  Lothe  (1968)  as 

ph  Xi(x^-xl) 

23c(1-v)  (^^  +  ^^2) 

Substituting  in  Eq.  (1)  by  the  magnitude  of  the  Burger's  vector  b,  it  is  seen  that  this  resolved 
shear  stress  caused  by  the  displacement  of  an  edge  dislocation  of  this  Burger's  vector  from  {-d,0) 
to  (d,0)  is  exactly  the  same  as  given  be  Eq.  (3).  Hence,  this  plastic  strain  in  this  thin  slice  can  be 
caused  by  the  displacement  of  dislocations.  This  can  also  be  interpreted  as  having  a  positive 
dislocation  at  (d,0)  and  a  negative  one  at  i-d,0). 

Consider  a  most  favorably  oriented  crystal  located  at  a  free  surface  of  an  aluminum 
polycrystal  loaded  in  tension.  All  crystals  are  taken  to  have  the  same  critical  shear  stress.  The 
crystal  with  the  highest  resolved  shear  stress,  known  as  the  most  favorably  oriented  crystal  will 
first  reach  the  critical  shear  stress  and  slides.  Under  uniaxial  tension,  this  crystal  has  a  slip 
direction  a  and  a  normal  p  to  the  slip  plane  making  45®  with  the  loading  axis.  A  small  uniform 
initial  resolved  shear  stress  is  assumed  to  exist  in  a  thin  slice,  which  will  slide  first.  This  slice 

is  divided  along  the  width  into  iV^-parallelogram  grids  Sn,  n  =  1,  2,  ...,  iV,  within  each  of  which 
this  plastic  strain  is  taken  to  be  uniform.  From  the  condition  that  the  resolved  shear  stress  Xap 

in  each  sliding  grid  equals  the  critical  shear  stress  x‘^,  the  slip  distribution  in  this  sliee  has  been 
determined  (Lin,  1972,  1992).  The  polycrystal  is  of  fine-grain  and  the  equivalent  force  may  be 
considered  as  acting  in  a  semi-infinite  solid.  The  thickness  of  the  slice  is  much  smaller  than  its 
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length.  The  plastic  strain  due  to  slip  is  assumed  to  be  constant  along  the  length,  so  the  deformation 
caused  by  the  equivalent  forces  are  taken  to  be  of  plane  strain.  The  plane  stress  solution  of  the 
stress  field  due  to  a  point  force  applied  in  a  semi-infinite  plate  has  been  given  by  Melan  (Tung  and 
Lin,  1966).  His  solution  was  modified  by  Tung  and  Lin  (1966)  for  plane  strain. 

Using  the  plane  strain  solution,  we  calculate  the  average  resolved  shear  stress  Xap  in  the 
grid  caused  by  a  unit  uniform  plastic  strain  in  the  n^  grid.  Denoting  this  average  stress 
by  Gmn^  we  can  express  the  residual  resolved  shear  stress  in  the  m^  grid  caused  by  plastic 

strain  in  all  grids  as 


n 


where  n  is  summed  over  all  grids  with  plastic  strain.  The  resolved  shear  stress  and  strain  in  this 
section  refer  to  the  ap-slip  system  only,  hence  the  subscripts  ap  are  dropped,  giving 

<  ■  2°-^^  <'*> 

n 


After  plastic  strain  occurs,  the  total  resolved  shear  stress  is  given  by 

i  r  a 
-t  =  X  +  X  +  X 

where  x“  is  the  resolved  shear  stress  caused  by  the  applied  load.  For  a  grid  to  start  or  continue 
sliding,  this  total  resolved  shear  stress  must  be  equal  to  the  critical  shear  stress 

X  =  x‘  +  tT  +  x“  =«  x*"  (5) 

Writing  this  in  incremental  form, 

Ax'^  +  Ax“  =  Ax*^  (6) 

Plcistic  deformation  in  metals  is  highly  locali2ed.  The  microscopic  plcistic  strain  in  the  thin 
slice  is  much  larger  than  the  macroscopic  plastic  strain  in  the  metal.  Hence,  the  strain  hardening  of 
the  thin  slice  is  much  less  than  that  of  the  metal  and  is  neglected.  This  gives  a  constant  critical 
resolved  shear  stress  Since  the  initial  shear  stress  does  not  change  with  loading,  with 
Ax*^  =  0 ,  we  have  in  any  sliding  grid 
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(7) 


=  0 


This  equation  is  ^plied  to  all  grids  with  incremental  plastic  strain.  We  have  as  many  equations  as 
the  unknowns.  Thus  the  plastic  strain  in  all  grids  were  solved  (Lin,  1972).  Multiple  slip  bands 
have  been  observed  in  the  most  favorably  oriented  crystal  at  free  surface.  For  numerical  study,  we 
consider  ten  slip  bands  spaced  at  3  pm  apart  and  with  initial  shear  stress  in  different  bands 
indicated  in  Figure  2.  The  polycrystal  is  taken  to  be  of  pure  aluminum  with  a  critical  shear  stress 
of  53.5  psi.  The  polycrystal  is  subject  to  a  tension  along  X2-direction.  The  calculated  plastic  strain 
distributions  of  the  different  bands  at  different  applied  loads  1:22  are  shown  in  Figure  3.  This  gives 
the  occurrence  of  shear  bands. 


Figure  2 

Most  favorably  oriented  crystal  with  ten  bands  and  initial  stress. 
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Surface  PLcstLc  Strain  In  Bands  xlQ 


Bands  {  Numbered  from  Bottom  to  Top  1 
Rgure  3 

Surface  plastic  strains  in  different  slip  bands. 


III.  Fatigue  Band  Formation 

More  than  90%  of  all  catastrophic  failures  of  machines,  vehicles,  and  other  structures 
occurring  in  practice  are  caused  by  fatigue  of  materials.  Fatigue  failures  occur  under  loadings 
substantially  lower  than  the  yield  strength  of  the  material.  Fatigue  cracks  occur  in  two  stages: 
crack  initiation  and  crack  propagation.  Crack  propagation  predominates  fatigue  life  in  low-cycle 
fatigue,  while  crack  initiation  predominates  the  life  in  high-cycle  fatigue  and  is  also  a  necessary 
stage  prior  to  crack  propagation.  Hence,  fatigue  crack  initiation  is  of  both  scientific  interest  and 
practical  importance. 

MacCammon  and  Rosenberg  (1957)  and  MacCrone  et  al  (1959)  showed  that  metals  are 
subject  to  failure  at  temperatures  as  low  as  1 .7®K.  This  indicates  that  although  surface  corrosion. 
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gas  adsorption,  gas  diffusion  into  the  metal,  and  vacancy  diffusion  to  form  voids  can  have  an 
important  effect  on  fatigue  deformation,  they  are  not  necessary  for  fatigue  failure.  This  seems  to 
leave  mechanics  (i.e.,  the  local  stress  and  strain),  as  the  basic  mechanism  of  fatigue. 

When  a  cyclic  loading  is  applied  to  a  metal,  the  individual  grains  begin  to  show  fine  slip 
markings  after  certain  cycles.  As  loading  continues,  some  of  the  lines  intensify  and  become  dark. 
These  intense  slip  bands  appear  to  be  the  source  of  fatigue  cracks.  If  the  surface  of  the  fatigued 
metal  is  electro-polished,  the  fainter  slip  markings  can  be  removed  but  not  the  darker  bands.  These 
are  known  cis  the  persistent  slip  bands.  The  life  of  a  specimen  can  be  prolonged  greatly  by 
removing  the  surface  layer  periodically.  A  microscopic  study  of  the  development  of  cracks  in 
copper  by  Thompson  et  al.  (1956)  verified  the  initiation  of  fatigue  cracks  in  these  bands. 


(a)  Intrusions  and  Extrusions 

Forsyth  and  Stubbington  (1955),  and  Hull  (1958)  found  thin  ribbons  of  metals  extruded 
from  die  slip  bands  during  fatigue  loading.  Extrusions  of  about  0.1  pm  thick  sometimes  reached  a 
height  of  10  pm  and  varied  in  width  from  about  1  pm  to  a  substantial  fraction  of  the  grain 
dimension.  This  formation  of  extrusions  clearly  are  associated  with  shears  of  opposite  signs  on 
the  opposite  sides  of  an  extrusion.  The  reverse  of  extrusion  (i.e.,  intrusion)  has  also  been 
observed.  Both  extrusions  and  intrusions  grow  monotonically  in  depth  and  width  with  cycles  of 
loading.  These  observations  are  pertinent  in  developing  the  theory  of  fatigue  crack  initiation. 

Fatigue  cracks  generally  initiate  at  a  free  surface.  We  consider  a  thin  slice  of  metal  in  a 
most  favorably  oriented  crystal  at  a  free  surface  of  a  polycrystal  subject  to  a  cyclic  tension  and 
compression.  The  slip  plane  and  slip  direction  of  this  crystal  form  45®  with  the  specimen  axis.  We 
consider  the  loading  to  be  of  low  amplitude.  Then  the  plastic  deformation  mainly  occurs  in  this 
most  favorably  oriented  crystal. 


(b)  A  Polycrystal  Fatigue  Crack  Initiation  Model 

Most  fatigue  cracks  initiate  at  a  free  surface.  To  relieve  the  same  amount  of  resolved  shear 
stress  in  a  thin  slice,  a  greater  amount  of  slip  (plastic  shear  strain)  is  required  near  the  free  surface 
than  at  the  interior  of  metals  (Lin  and  Ito,  1967).  Hence,  the  thin  slice  subject  to  alternate  forward 
and  reversed  loadings  of  this  model  is  taken  to  be  in  the  most  favorably  oriented  crystal  at  a  free 
surface  of  a  polycrystal,  as  shown  in  Figure  4.  This  polycrystal  is  loaded  in  cyclic  tension  and 
compression  of  low  amplitude.  Plastic  deformation  occurring  in  this  most  favorably  oriented 
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crystal  in  the  polycrystal  causes  slip  to  occur  in  the  most  favorable  slip  system  to  form  a  persistent 
slip  band. 


Figure  4 

Most  favorably  oriented  crystal  at  free  surface. 


The  thickness  of  the  slices  P,  Q,  suid  R  is  much  smaller  than  the  length  (dimension  along 
xs-axis).  The  plastic  strain  caused  by  slip  is  assumed  to  be  constant  along  this  length.  Hence,  the 
equivalent  forces  caused  by  this  plastic  strain  are  considered  to  act  in  a  semi-infinite  elastic  medium 
under  plane  deformation. 


(c)  Initial  Resolved  Shear  Stress 

Initial  defects  always  exist  in  metals  and  cause  an  initial  stress  field  x* .  During  loading, 
when  the  resolved  shear  stress  in  some  region  reaches  the  critical  shear  stress  x'",  slip  occurs. 
After  unloading,  this  slip  remains  and  induces  a  residual  resolved  shear  stress  x\  The  resolved 
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shear  stress  caused  by  loading  is  denoted  by  .  Hence,  the  total  resolved  shear  stress  is  the  sum 
of  these  three  stresses. 


T  =  x‘  +  /  +  x“  (8) 

For  an  extrusion  to  initiate,  positive  shear  has  to  occur  in  a  thin  slice  P  and  negative  shear  in  a 
closely  located  slice  Q  (see  Figure  4).  The  initial  shear  stress  field  x'  favorable  to  this  sequence  of 
slip,  clearly  is  one  having  positive  shear  stress  in  P  and  negative  in  Q.  Such  an  initial  stress  field 
can  be  provided  by  an  initial  tensile  strain  in  R.  (The  repetition  of  Greek  subscripts  does  not 
denote  summation.)  This  positive  can  be  provided  by  a  row  of  interstitial  dislocation  dipoles 
and  a  negative  by  vacancy  dipoles  as  suggested  by  Lin  and  Ito  (1969a).  Recently, 
Antonopoulos  etal.  (1967)  and  Mughrabi  et  al.  (1983)  have  shown  that  the  ladder  structure  in  a 
persistent  slip  band  (PSB)  can  be  represented  by  an  array  of  dislocation  dipoles,  causing  initial 
resolved  shear  stresses  at  the  interface  between  the  PSB  and  the  matrix. 


(d)  Gating  Mechanism  Provided  by  Microstresses 

A  tensile  loading  X22  on  the  polycrystal  (Figure  4)  produces  a  positive  x“  in  the  whole 
crystal.  Taking  x‘  to  be  positive  in  P  and  negative  in  Q,  we  have  x'  +  x"  in  P  reaching  the  critical 
shear  stress  first,  and  hence  P  slides.  The  plastic  strain  caused  by  slip  is  taken  to  be  constant 
along  JC3-axis,  hence  dx^/dx^  =  0.  The  equilibrium  condition  with  no  body  force,  gives 


^gg  ,  _  0 

da  dp 


(9) 


Since  dx^Jda  is  finite,  dx^fd^  must  be  also  finite,  and  the  changes  in  Xap  across  the  small 

distance  between  P  and  Q  is  very  small.  Therefore,  the  slip  in  P  relieves  not  only  the  positive 
shear  stress  in  P,  but  also  in  its  neighboring  region,  including  Q.  The  relief  of  positive  shear  stress 
is  the  same  as  the  increase  of  negative  resolved  shear  stress  in  Q  to  cause  Q  to  slide  more  readily  in 
the  reverse  loading.  The  negative  slip  in  Q  during  the  reverse  loading  relieves  the  negative  shear 
stress  not  only  in  Q,  but  also  in  P,  thus  causing  P  more  readily  to  slide  in  the  next  forward  loading. 

This  process  is  repeated  for  every  cycle  thus  providing  a  gating  mechanism  for  a  monotonic 
buildup  of  local  slip  strain  in  P  and  Q,  pushing  R  out  of  the  free  surface  and  starts  an 

extmsion.  Interchanging  the  signs  of  the  initial  stresses  in  P  and  Q  initiates  an  intrusion  instead  of 
an  extrusion.  This  theory  is  extensively  supported  by  metallurgical  observations  (Lin,  1977). 
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This  gating  mechanism  was  introduced  by  Lin  and  Ito  (1969b)  and  was  given  in  dislocation  arrays 
by  Tanaka  and  Mura  (1981)  and  is  also  discussed  in  Suresh's  Fatigue  of  Materials  (1991). 


(e)  Activation  of  Second  Slip  System 

The  buildup  of  this  local  slip  strain  e'^  in  P  and  Q  is  caused  by  the  positive  and  negative 
initial  shear  stress  which,  in  turn,  is  caused  by  in  R.  If  R  were  cut  out,  this  free  length  of 

R  would  be  longer  than  the  slot  cut,  by  an  amount  referred  to  as  the  "static  extrusion" 
(Antonopoulos  etal.,  1967;  Mughiabi  etal.,  1983).  This  e‘^  causes  an  initial  compression  in 

R.  As  the  extrusion  grows  under  cyclic  loading,  the  slice  R  increases  in  length.  This  elongation 
causes  the  compression  to  decrease.  A  question  has  been  raised  as  to  whether  the  extrusion 
growth  will  ce2ise  after  extrusion  has  reached  the  static  extrusion.  For  aluminum  and  its  alloys,  the 
residual  tensile  stress  caused  by  elongation  in  R  due  to  extrusion  can  cause  changes  of 

resolved  shear  stresses  in  all  12  slip  systems  of  the  face-centered  cubic  crystal.  The  resolved  shear 
stress  in  one  slip  system  may  reach  the  critical  shear  stress  and  slide.  The  plastic  strain  caused 

by  slip  in  this  system  has  a  tensor  component  e^  just  like  e^^^  in  causing  the  positive  and  negative 
,  respectively,  in  P  and  Q.  This  secondary  slip  has  been  shown  by  Lin  et  al.  (1989)  to  increase 
greatly  the  extent  of  extrusion  and  intrusion.  This  causes  equivalent  forces  in  addition  to 
Fj  and  F2.  With  F3,  the  strains  in  the  slices  are  no  longer  of  plane  strain  and  hence  are  taken  to 
be  of  generalized  plane  strain;  i.e.. 


Ui=Ui(Xi,X2)  1,2,3 


(10) 


The  solution  of  a  similar  problem  was  shown  in  the  analysis  of  prismatic  anisotropic  bars  by 
Lekhnitski  (1963).  The  stress-strain  relation  gives 


Cj-  =  2ix 


1  -  V 


•6«e  +  - 


1 1  du,  du,^ 


2\dx, 


dXi, 


(11) 


where 


^  diu  5m, 
9  =  — i-  +  — - 

dxi  6x2 
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and  V  is  the  Poisson's  ratio, 
yields 


Substitution  of  the  expression  into  the  conditions  of  equilibrium 


v\.  + 


1  ae 

+ 

1  -  2v  dXi 


i  =  1,2 


(12) 


and 


=  0 


(13) 


where 


dx^  ^  dx^ 


and  Fi  and  F3  can  be  either  the  applied  forces  or  the  equivalent  forces. 

Equations  ( 12)  and  ( 13)  are  not  coupled  and  can  be  solved  separately.  Lin  and  Lin  ( 1983) 
have  shown  the  solution  of  this  stress  field  caused  by  F,  in  Eq.  (12)  by  using  Airy's  stress 
functions.  An  initial  shear  stress  field  in  three  thin  slices  P,  Q,  and  R  in  Figure  1  caused  by 

varying  linearly  from  zero  at  the  free  surface  to  a  maximum  at  the  interior  boundary  in  the  outer 
two  slices,  P  and  Q,  and  by  an  of  opposite  sign  and  of  double  magnitude  in  the  center  slice  R 

was  numerically  calculated  by  Lin  and  Lin  (1983).  These  calculated  initial  shear  stresses  were 
found  to  be  uniform  in  the  two  outer  slices,  one  is  positive  and  the  other  is  negative.  The  regions 
outside  of  these  two  slices  have  a  negligible  small  amount  of  resolve  shear  stresses.  This  initial 
stress  field  so  derived  from  an  initial  strain  clearly  is  favorable  to  the  initiation  of  extrusion  and 
satisfies  both  the  equilibrium  and  compatibility  conditions.  This  shows  that  initial  resolved  shear 
stresses  in  P  and  Q  can  be  provided  by  an  initial  strain  and  satisfy  equilibrium  conditions  (Lin 
and  Lin,  1983).  To  solve  Eq.  (13)  for  F3,  we  write  ^>3(0:,  J)  as  pM3(x,J).  Ui(x,x)  represents 
the  displacement  «3  at  points  due  to  force  applied  at  x.  The  stress  is  then 

tijCx.T)  =  ,  'C23(^>^)  =  (14) 

For  a  unit  concentrated  force  F3  at  T ,  Eq.  ( 13)  gives 

V^<I>3 +6(x-3t)-0  (15) 

where  6(Ar -  x)  is  the  Dirac  Delta  function.  With  the  boundary  conditions  of 
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00 


5^0  I  2  ^ 

— ^  =  0  at  JCj  =  0  and  — ^  -»  0  as  +X2 

dx^  8X2 

we  have 


4>3(x,3c)  =  -— (inXi  +  InXj  (16) 

4ji 

where 

Xi  +(^X2-X2f 

Xj  =(Afi  +  ^)^  +{X2-^f 


Hence  the  stress  fields  caused  by  the  equivalent  forces  due  to  slip  strain  can  be  readily  be 
calculated.  This  micromechanic  theory  of  high-cycle  fatigue  crack  initiation  has  extensive 
metallurgical  verifications  (Lin,  1992). 


IV.  Interaction  of  Multiple  Fatigue  Bands 

In  the  early  stages  of  high-cycle  fatigue,  the  fatigue  bands  are  widely  spaced  as  light  slip 
markings  (Kennedy,  1963)  in  the  most  favorably  oriented  crystal  at  a  free  surface  of  a  polycrystal. 
The  interaction  of  these  bands  is  small,  so  the  slip  distributions  in  the  band  can  be  analyzed  as  a 
single  fatigue  band  in  the  crystal.  As  the  fatigue  loading  proceeds,  more  slip  lines  appear  and  the 
interaction  of  neighboring  fatigue  bands  in  the  crystal  needs  to  be  considered.  Since  the  crystals 
surrounding  this  most  favorable  one  are  less  favorably  oriented,  these  multiple  fatigue  bands  can 
occur  in  this  most  favorable  crystal,  while  no  fatigue  band  in  the  surrounding  crystals. 

One  case  of  the  interaction  of  slip  bands  in  HCF  crack  initiation  was  analyzed  by  Lin  and 
Chen  (1989).  These  bands  were  assumed  to  exist  in  a  most  favorably  oriented  crystal  at  a  free 
surface  of  a  polycrystal.  Each  band  is  taken  to  consist  of  3  thin  slices  P,  Q,  and  R  as  shown  in 
Figure  4.  P  and  Q  in  each  band  are  assumed  to  have  equal  and  opposite  initial  resolved  shear 
stresses.  These  initial  stresses  are  different  in  different  bands.  Plastic  strain  distributions  in  these 
bands  under  cyclic  loadings  were  calculated.  The  maximum  shear  strain  at  the  free  surface  is  an 
indication  of  extrusion  height  and/or  intrusion  depth,  which  is  taken  as  a  measure  of  crack 
initiation.  This  indicates  the  level  of  fatigue  damage.  The  applied  cyclic  sheeu  stress  versus  the 
number  of  cycles  to  yield  different  amount  of  surface  shear  strains  were  calculated.  The  shape  of 
these  curves  is  similar  to  that  of  Wohler's  diagrams  and  that  of  Coffin-Mansan  equation.  The 
results  seem  to  give  good  physical  explanations  to  these  diagrams  and  equations. 
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V.  Propagation  of  Fatigue  Bands  across  Grain  Boundaries 


Grain  boundaries  are  only  few  atomic  spacings.  These  boundaries  are  assumed  to  have 
zero  thickness.  Fatigue  fracture  of  a  polycrystal  is  either  transgranular  or  intergranular.  This  mode 
of  fracture  affects  the  ductility  of  the  metal.  It  is  expected  that  the  propagation  of  fatigue  band 
across  the  grain  boundary  more  likely  yields  transgranular  crack  initiation  and  the  stopping  of  the 
propagation  at  the  grain  boundary  tends  to  cause  an  intergranular  crack.  A  single  fatigue  band  in 
the  most  favorably  oriented  crystal  at  a  free  surface  of  a  polycrystal  has  been  analyzed.  The 
dimensions  of  this  fatigue  band  in  surface  crystal  and  its  neighboring  crystal  is  shown  in  Figure  5. 
This  polycrystal  is  a  pure  aluminum  with  a  critical  shear  stress  of  53.3  psi.  An  initial  shear  stress 
x'  of  0.5  psi  in  P  and  Q  was  assumed.  The  polycrystal  was  subjected  to  a  cyclic  tension  and 
compression  of  107.7  psi. 


Figure  5 

Fatigue  band  across  the  grain  boundary. 


The  slip  directions  and  slip  planes  of  a  f.c.c.  crystal  are  shown  in  Figure  6.  Let  a2 
correspond  to  the  ap  slip  system  in  Figure  4,  which  is  referred  to  as  the  primary  slip  system.  The 
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extrusion  growth  causes  the  tensile  stressxaa  in  R.  There  are  four  slip  system  ci,  C3,  di,  and  d2 
equally  favorable  under  Xaa-  Of  these  four,  C3  is  most  favorable  under  X22,  and  hence  will  be 
activated  next  and  is  referred  to  as  the  second  active  slip  system  (Lin,  1992).  The  calculated  plastic 
strain  distribution  of  the  primary  slip  system  versus  the  distance  from  the  free  surface  at  different 
cycles  of  loading  are  shown  in  Figure  7.  The  plastic  strain  distributions  of  the  second  slip  system, 
C3,  in  R  are  shown  in  Figure  8.  The  surface  plastic  strain  approaches  a  saturated  value  at  about 
650  cycles.  It  is  seen  that  very  little  second  slip  in  the  neighboring  crystal.  This  propagation 
mainly  depends  on  the  difference  of  orientation  of  the  two  crystals. 


Crystallographic  directions  of  a  f.c.c.  crystal. 


Multiple  fatigue  bands  have  been  observed  on  the  surface  crystal.  As  the  magnitude  of  the 
alternate  loading  increases,  those  crystals  surrounding  the  surface  crystal  will  also  plastically 
deform  and  give  residual  stress  fields,  which  will  interact  with  fatigue  bands  of  the  surface  crystal. 
The  resultant  slip  distributions  in  these  bands  have  been  analyzed  by  Lin  and  Chen  (1989). 
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Rgure  7 

Plastic  strain  distribution  of  the  primary  slip  system  in  P. 


'c>  a 


Rgure  8 

Plastic  strain  distribution  of  the  second  slip  system  in  R. 
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VI.  HCF  Crack  Initiation  of  Intermetallic  Compounds 


Ordered  intermetallic  compounds  favor  planar  slip  eind  tend  to  have  no  cross-slip  or 
multiple  slip.  This  is  due  to  the  long  Burger's  vector  of  these  intermetallics.  Test  results  show  that 
long-rang  order  in  intermetallic  compounds  substantially  increase  their  fatigue  lives  in  stress- 
controlled  cycling  but  decreases  in  the  strain-controlled  fatigue  resistance.  Stoloff  et  al.  (1986) 
indicated  that  the  least  understood  aspects  of  the  cyclic  behavior  of  intermetallic  compounds  are  the 
processes  leading  to  crack  initiation.  Using  the  inhibition  of  cross  or  multiple  slip  in  the  ordered 
compounds,  the  present  micromechanic  model  has  been  found  to  satisfactorily  explain  the  cyclic 
behavior  of  the  intermetallic  compounds  (Lin  and  Chen,  1992, 1993). 


VII.  Effect  of  Crystal  Elastic  Anisotropy  on  Polycrystal 

Fatigue  Band 

Previous  studies  on  fatigue  bands  were  for  aluminum  and  its  alloys.  The  anisotropy  of 
their  elastic  constants  is  small  and  accordingly  neglected.  However,  this  anisotropy  of  other  metals 
such  as  titanium  and  some  intermetallic  alloys  is  not  small.  Hence,  in  this  section,  anisotropy  is 
here  analyzed. 

In  an  earlier  study,  the  anisotropy  of  the  most  favorably  oriented  crystal  at  a  free  surface  is 
considered.  The  metal  exterior  to  the  surface  is  assumed  to  be  elastically  isotropic  and 
homogeneous.  Thus,  it  is  reduced  to  an  inclusion  problem.  This  heterogeneous  solid  is 
transformed  into  a  homogeneous  one  by  using  Eshelby's  Equivalent  Inclusion  Method  (Eshelby, 
1957,  1961,  Mura,  1982).  In  this  method,  an  extra  set  of  eigenstrains,  which  behave  the  same  as 
the  inelastic  strains,  is  introduced.  The  eigenstrains  explicitly  depend  on  the  applied  stresses  and 
plastic  strains  and  thus  can  be  implicitly  incorporated  in  the  formulation  of  solutions.  The  surface 
crystal  is  assumed  to  have  the  anisotropic  elastic  constants  of  NisAl  monocrystalline  intermetallic 
compound  (Y ang,  1985).  The  surrounding  crystals  are  assumed  to  have  isotropic  elastic 
constants.  The  plastic  strain  distributions  and  the  cumulative  surface  plastic  strain  in  the  fatigue 
band  versus  the  number  of  loading  cycles  are  calculated,  and  the  effect  of  elastic  anisotropy  on  the 
formation  of  extrusions  and  intrusions  is  shown.  The  details  of  the  analysis  are  shown  in  the 
paper  published  in  AMD-Vol.  204,  Numerical  Methods  in  Structural  Mechanics,  ASME,  1995.  A 
copy  of  this  paper  is  given  in  the  appendix. 
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In  the  above  study,  the  crystals  surrounding  the  surface  crystal  were  assumed  to  be 
isotropic.  In  a  later  refined  study,  the  orientations  of  the  surrounding  crystal  were  considered. 
This  introduces  different  anisotropic  elastic  constants  to  different  crystal  referring  to  the  specimen 
axes.  This  model  is  more  realistic.  The  plastic  shear  strain  distribution  in  the  fatigue  bands  were 
calculated.  This  cumulative  surface  plastic  shear  strains  versus  the  number  of  loading  cycles  for 
the  case  with  the  consideration  of  crystal  anisotropy  and  for  the  case  without,  were  calculated  and 
shown.  It  is  shown  that  the  neglect  of  this  crystal  anisotropy  may  cause  significant  error  in 
calculating  the  slip  distribution  and  in  estimating  the  fatigue  crack  initiation.  The  analysis  is 
published  in  the  Journal  of  Engineering  Materials  and  Technology  (Teng  and  Lin,  1995).  A  copy 
of  the  paper  is  given  in  the  appendix. 
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Abstract — Tests  have  shown  that  ordered  intermetallic  compounds  favor  planar  slip  and  tend  to  have 
no  cross-slip  or  multi-slip.  This  was  explained  by  the  long  Burger’s  vector  in  these  ordered  compounds. 
The  flow  stress  of  the  ordered  has  been  found  to  be  less  than  that  of  the  disordered  ones.  However  the 
high-cycle  fatigue  life  of  the  ordered  was  found  to  be  significantly  more  than  that  of  the  disordered.  This 
is  believed  to  be  due  to  low  rate  of  crack  initiation.  Cracks  have  been  found  to  initiate  exclusively  at  slip 
bands  extrusions  and  intrusions  in  intermetallic  compounds.  This  paper  shows  a  micromechanic  fatigue 
model.  This  model  was  used  to  calculate  the  extrusion  growth  of  a  single  fatigue  band  and  then  multiple 
fatigue  bands  in  the  most  favorably  oriented  crystal.  It  was  found  that  extrusion  growth  in  the  ordered 
alloy  is  significantly  less  than  that  in  the  disordered  ones.  The  curves  to  give  given  plastic  shear  strains 
of  the  extrusion  at  the  free  surface  were  calculated  for  both  the  ordered  and  disordered  alloys.  It  was  found 
that  the  ordered  alloy  has  significantly  higher  S-N  curves  as  observed  in  experiments. 

Resume — Des  essais  monlrent  que  les  composes  intermetalliques  ordonnes  favorisent  un  glissement  plan 
et  tendent  a  ne  pas  operer  de  glissements  devies  simple  ou  multiple.  Ceci  s’explique  par  le  long  vecteur 
de  Burgers  des  ces  composes  ordonnes.  On  trouve  que  la  contrainte  d’ecoulement  plastique  des  materiaux 
ordonnes  est  inferieure  a  celle  des  materiaux  desordonnes.  Cependent,  on  trouve  que  la  duree  de  vie  en 
cyclage  a  frequence  elevee  des  materiaux  ordonnes  est  beaucoup  plus  grande  que  celle  des  materiaux 
deordonnes.,On  pense  que  ceci  est  du  a  la  faible  vitesse  d’initiation  des  fissures.  On  trouve,  dans  les 
composes  intermetalliques,  que  les  fissures  commencent  exclusivement  aux  extrusions  et  intrusions  des 
bandes  de  glissement.  Cet  article  donne  un  modele  de  fatigue  micromecanique.  Ce  modele  est  utilise  pour 
calculer  la  croissance  de  Textrusion  d’une  bande  unique  de  fatigue,  puis  de  bandes  multiples  dans  un  cristal 
oriente  de  la  maniere  la  plus  favorable.  On  trouve  que  la  croissance  d’une  extrusion  dans  un  alliage 
ordonne  est  beaucoup  moins  importante  que  dans  des  alliages  desordonnes.  Les  courbes  S-Ny  destinees 
a  donner  des  deformations  plastiques  donnees  par  cisaillement  de  I’extrusion  sur  la  surface  libre  sont 
calculees  a  la  fois  pour  des  alliages  ordonne  et  desordonne.  On  trouve  que  Talliage  ordonne  a  des  courbes 
S-N  beaucoup  plus  hautes  que  celles  que  Ton  observe  experimentalement. 

Zusammenfassung — Versuche  haben  gezeigt,  daB  geordnele  intermetallische  Verbindungen  ebene  Gleitung 
bevorzugen,  also  Quergleitung  und  Vielfachgleitung  verraeiden.  Dieses  Verhalten  wird  iriit  dem  langen 
Burgersvektor  erklart.  Die  FlieBspannung  der  geordneten  Legierungen  ist  kleiner  als  die  der  entordneten. 
Allerdings  ist  die  Ermudungslebensdauer  der  geordneten  Legierungen  betrachtlich  groBer  als  die  der 
entordneten.  Das  wirdeiner  niedrigen  Ri Child ungsrate  zugeschrieben.  Die  Risse  entstehen  ausschlieClich 
an  Gleitbandextrusionen  und  -intrusionen.  Diese  Arbeit  baut  ein  mikromechanisches  Ermudungsmodell 
auf.  Mit  diesem  Modell  wird  das  Extrusionswachstum  eines  einzelnen  Ermiidungsbandes  und  dann  vieler 
Ermudungsbander  in  den  am  gunstigsten  orientierten  Kornern  berechnet.  Es  ergibt  sich,  daC  das 
Extrusionswachstum  in  den  geordneten  Legierungen  bedeutend  gerinfeer  ist  als  in  den  entordneten.  Die 
5-A^-Kurven,  welche  die  plastischen  Scherungen  der  Extrusionen  an  der  freien  Oberflache  angeben, 
werden  fiir  geordnele  und  ungeordnete  Legierungen  berechnet.  Das  Ergebnis  ist,  daC  die  geordnete 
Legierung  betrachtlich  hohere  S-A^-Kurven  aufweist  als  gemessen  im  Experiment. 


1.  INTRODUCTION 

The  ordered  structure  of  intermetallic  compounds 
gives  a  dislocation  Burgers  vector  quite  long  com¬ 
pared  to  that  of  disordered  alloys.  The  dislocation 
with  this  long  Burger’s  vector  is  known  as  a  super¬ 
dislocation.  In  a  f.c.c.  material,  the  Burger’s  vector  is 
j  [1 10],  while  in  an  ordered  intermetallic  compounds 
such  as  NijAl,  the  Burger’s  vector  is  [110].  The 
increased  length  of  the  Burger’s  vector  in  inter¬ 
metallic  compounds  favors  planar  slip  and  tends  to 
inhibit  cross-slip  and/or  multiple  slip. 


Fatigue  tests  of  FeCo-V  show  that  slip  along  (1 10) 
occurs  for  both  ordered  and  disordered  alloy 
samples.  Cracks  initiated  almost  exclusively  at  slip 
bands,  intrusions  and  extrusions.  Tests  carried  out 
in  high  vacuum  (5  x  torr)  also  have  shown  that 
slip  bands  are  preferred  sites  for  crack  initiation 
and  propagation  for  both  ordered  and  disordered 
conditions.  Similar  tests  on  Ni3Al  show  the  develop¬ 
ment  of  extrusions  and  intrusions  into  cracks. 
Hence  the  height  of  the  extrusion  or  the  depth  of 
the  intrusion  is  here  taken  as  a  measure  of  fatigue 
damage. 
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This  influence  of  long  range  order  on  stress- 
controlled  high  cycle  fatigue  tested  at  25’"C  was  first 
reported  for  Ni^Mn  and  FeCo-V.  In  both  cases 
high-cycle  fatigue  lives  were  substantially  increased, 
in  spite  of  a  decrease  in  yield  stress  with  order  for 
FeCo-V  [2, 3].  Later,  high  cycle  fatigue  tests  were 
extended  to  FeNi^-V  polycrystals,  NijAl  single 
crystals  etc.  [4].  In  both  cases,  fatigue  lives  were 
substantially  increased  in  spite  of  a  decrease  in  flow 
stress  with  order  for  FeCo-V.  It  was  found  that  the 
ratio  of  the  fatigue  limit  to  yield  strength  exceeds  0.5 
for  most  of  the  ordered  alloys,  while  in  the  disordered 
condition,  the  ratio  decreases  significantly,  [5].  Since 
crack  propagation  is  much  more  rapid  in  ordered 
conditions  as  indicated  by  Stolofif  et  al  [5]  it  has  been 
concluded  that  ordering  substantially  delays  crack 
initiation  while  not  changing  the  location  of  the 
earliest  crack  [5].  They  further  indicated  that  the  least 
understood  aspects  of  cyclic  behavior  of  intermetallic 
compounds  are  the  processes  leading  to  crack 
initiation.  The  present  paper  is  an  attempt  to  give  the 
mechanism  of  high-cycle  fatigue  crack  initiation  of 
some  long  range  ordered  intermetallic  compounds. 
The  present  study  analyzes  the  micro  stress  and  strain 
fields  of  the  fatigue  bands.  These  stress  and  strain 
fields  show  analytically  the  long  range  order  effect 
on  fatigue  crack  initiation.  This  analysis  is  not  for 
low  cycle  fatigue.  As  indicated  by  Stoloff  et  al  [5], 
there  is  no  evidence  that  ordering  is  beneficial  under 
strain-controlled  low  cycle  fatigue  conditions. 

Flow  stress  of  many  intermetallic  alloys  increase 
with  quenched  temperature  to  a  peak  temperature  [2]. 
Octahedral  slip  was  found  to  be  predominant  at 
temperatures  below  the  peak  temperature,  while 
primary  cube  slip  is  prevalent  above  the  peak  tem¬ 
perature  [4].  The  present  calculated  results  are  to  be 
compared  to  the  test  results  conducted  at  25°C,  hence 
octahedral  slip  is  considered.  In  the  present  analysis, 
the  intermetallic  compounds  are  assumed  to  follow 
Schmids  law. 

2.  A  POLYCRYSTAL  FATIGUE  MODEL 

Extrusions  and  intrusions  in  fatigue  specimens 
were  first  discovered  by  Forsyth  and  Stubbington  [6] 
and  later  found  by  a  number  of  other  physicists  and 
material  scientists.  These  extrusions  and  intrusions 
are  favourable  sites  for  fatigue  cracks.  Following 
the  clue  provided  by  the  observation  on  extrusion 
and  intrusions,  a  number  of  theories  of  fatigue  crack 
initiation  were  proposed  by  different  investigators  [7]. 
One  theory  considered  a  column  of  metal  containing 
a  single  screw  dislocation  traveling  a  complete  circuit 
and  the  volume  contained  in  the  circuit  is  translated 
parallel  to  the  dislocation,  this  causes  the  metal  to 
extrude.  This  mechanism  does  not  explain  why  the 
dislocation  under  cyclic  stressing,  does  not  oscillate 
back  and  forth  along  the  same  path  rather  than 
traversing  a  closed  circuit.  Drawbacks  of  this  and 
other  theories  are  discussed  by  Kennedy  [7].  Clearly, 


\  \  \  \  \  \ 
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Fig.  1.  Most  favorably  oriented  crystal  at  free  surface. 

some  gating  mechanism  is  needed  to  convert  the 
back  and  forth  oscillation  of  the  dislocation  into  a 
unidirectional  circuit. 

Most  fatigue  cracks  initiate  at  a  free  surface.  To 
relieve  the  same  amount  of  resolved  shear  stress  in  a 
thin  slice,  a  greater  amount  of  slip  (plastic  shear  strain) 
is  required  near  the  free  surface  than  at  the  interior 
of  metals  [8].  Hence,  the  thin  slice  subject  to  alternate 
forward  and  reversed  loadings  of  this  model  is  taken 
to  be  in  the  most  favorably  oriented  crystal  at  a  free 
surface  of  a  polycrystal,  as  shown  in  Fig.  1.  This  poly¬ 
crystal  is  loaded  in  cyclic  tension  and  compression 
of  low  amplitude  and  plastic  deformation  essentially 
occurs  only  in  the  most  favourably  oriented  crystal. 
Figure  1  corresponds  to  a  persistent  slip  band, 
which  was  formed  after  slip  had  occurred  in  the  most 
favourable  slip  system  in  the  whole  crystal. 

(a)  Initial  resolved  shear  stress 

Initial  defects  always  exist  in  metals  and  cause 
an  initial  stress  field  t\  During  loading,  when  the 
resolved  shear  stress  in  some  region  reaches  the 
critical  shear  stress  slip  occurs.  After  unloading, 
this  slip  remains  and  induces  a  residual  resolved  shear 
stress  tL  The  resolved  shear  stress  caused  by  loading 
is  denoted  by  t“.  Hence  the  total  resolved  shear  stress 
is  the  sum  of  these  three  stresses 

T  =  T* -f- -h  T®.  (1) 

For  an  extrusion  to  initiate,  positive  shear  has  to 
occur  in  a  thin  slice  P  and  negative  shear  in  a  closely 
located  slice  Q  (Fig.  1).  The  initial  shear  stress  field 
t‘  favorable  to  this  sequence  of  slip,  clearly  is  one 
having  positive  shear  stress  in  P  and  negative  in  Q. 
Such  an  initial  stress  field  can  be  provided  by  an 
initial  tensile  strain  in  R,  (The  repetition  of 
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Greek  superscripts  does  not  denote  summation. )  This 
positive  can  be  provided  by  a  row  of  interstitial 
dislocation  dipoles  and  a  negative  e!,,  by  vacancy 
dipoles  as  suggested  by  Lin  and  Ito  [10].  Recently, 
Antonopoulos  et  aL  [11]  and  Mughrabi  et  ai  [12] 
have  shown  that  the  ladder  structure  in  a  persistent 
slip  band  (PSB)  can  be  represented  by  an  array  of 
dislocation  dipoles,  causing  initial  resolved  shear 
stresses  at  the  interface  between  the  PSB  and  the 
matrix. 


(b)  Gating  mechanism 


A  tensile  loading  X22  the  polycrystal  (Fig.  1) 
produces  a  positive  in  the  whole  crystal.  Taking 
t‘  to  be  positive  in  P  and  negative  in  Q,  we  have 
T*  -f  T*  in  jP  reaching  the  critical  shear  stress  r®  first; 
and  hence,  P  slides.  The  plastic  strain  caused  by  slip 
is  taken  to  be  constant  along  X3-axis,  hence 
dx^  =  0.  The  equilibrium  condition  with  no  body  force, 
gives 


da 


=  0. 


Since  dx^Jda  is  finite,  dx^^ldfi  must  be  also  finite,  and 
the  change  in  across  the  small  distance  between  P 
and  Q  is  very  small.  Therefore  the  slip  in  P  relieves 
not  only  the  positive  shear  stress  in  but  also  in  its 
neighboring  region  including  Q.  The  relief  of  positive 
shear  stress  is  the  same  as  the  increase  of  negative 
shear  stress.  Hence  this  slip  increases  the  negative 
resolved  shear  stress  in  Q  to  cause  Q  to  slide  more 
readily  in  the  reverse  loading.  The  negative  slip  in  Q 
relieves  the  negative  shear  stress  not  only  in  2,  but 
also  in  P,  thus  causing  P  more  readily  to  slide  in  the 
next  forward  loading.  This  process  is  repeated  for 
every  cycle  thus  providing  a  natural  gating  mechan¬ 
ism  for  a  monotonic  buildup  of  local  slip  strain  in 
P  and  pushing  R  out  of  the  free  surface  and  starts 
an  extrusion.  Interchanging  the  signs  of  the  initial 
stresses  in  P  and  Q  initiates  an  intrusion  instead  of 
an  extrusion.  This  theory  is  extensively  supported  by 
metallurgical  observations  [8]. 


3,  CALCULATION  OF 
RESIDUAL  STRESS  HELD 

The  resolved  shear  stress  is  the  sum  of  the  initial, 
applied  and  residual  shear  stresses  [equation  (1)]. 
To  calculate  the  residual  stress,  the  analogy  between 
plastic  strain  and  applied  force  [13]  is  applied.  It  has 
been  shown  that  the  equivalent  body  force  per  unit 
volume  along  the  X/-axis  due  to  plastic  strain  e*!j  in  a 
body  with  isotropic  elastic  constants  is 

(2) 

where  k  and  G  are  Lame’s  constants.  The  repetition 
of  an  italic  subscript  denotes  summation  and  the 
subscript  after  a  comma  denotes  differentiation,  with 
respect  to  the  coordinate  variable.  The  equivalent 


surface  force  per  unit  area  along  the  x,-axis  has  been 
shown  as 

=  (3) 

where  Vy  is  the  cosine  of  the  angle  between  the  exterior 
normal  to  the  surface  and  the  x,-axis. 

From  a  given  cj,  we  obtain  iv  and  S,  from 
equations  (2)  and  (3).  Then  we  apply  this  /;  and  5,  to 
the  body  and  obtain  the  strain  distribution  in  this 
body  by  elastic  solution.  The  residual  stress  field  in 
this  body  is  then 

Tly  =  (2  +  2Ge,y  -  X  ~  2Ge".)  (4) 

where  is  the  Kronecker  delta.  This  analogy  reduces 
to  Duhamel’s  analogy  for  thermal  stress  [12],  if  the 
plastic  strain  is  replaced  by  thermal  strain  ej. 

The  lengths  of  slices  P,  Q  and  R  along  the  X3- 
direction  are  much  larger  than  the  thickness  and  their 
inclined  lengths  as  shown  in  Fig.  1.  Hence,  strains  in 
the  major  central  portion  of  the  slices  are  taken  to  be 
of  generalized  plane  strain,  i.e. 

=  f  =  l,2,3.  (5) 

The  solution  of  a  similar  problem  was  shown  in  the 
analysis  of  prismatic  anisotropic  bars  by  Lekhnitski 
[14].  The  stress-strain  relation  gives 

^.7  =  +  \  (Uij  +  M;.,) j  (6) 

where  6  =  Wj ,  -h  1/2.2  and  v  is  Poisson’s  ratio. 
Substitution  of  the  expression  into  the  conditions 
of  equilibrium  yields 


and 


1  d0  Ff 
\-2vdx,'^  G 


=  0, 


/  =  1,2 


where 


dx\ 


(7) 

(8) 


Fi  and  F3  can  be  either  the  applied  forces  or  the 
equivalent  forces. 

Equations  (7)  and  (8)  are  not  coupled  and  can  be 
solved  separately.  Lin  [13]  have  shown  the  solution 
of  this  stress  field  caused  by  Py  in  equation  (7),  by 
using  Airy’s  stress  functions.  Let  t^(x,  x)  be  the  stress 
at  X  due  to  a  unit  force  along  x^-axis  applied  at  x 
of  this  semi-infinite  medium  and  <f>k0^,x)  be  the 
corresponding  Airy  stress  function,  where  x  denotes 
(xj,  X2,  X3).  Lin  and  Lin  [15]  have  expressed  the  stress 
components  in  terms  of  the  stress  function  as 


Tt,(x,S)  = 


dxl  ’ 


T^(X,X) 


T^(X,  X)  = 


6x2  ’ 


k  =  l,2 


(9) 
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|"(A,(x.x)=  -ip  -\-qKx:,~x^){0,  + 62) 
+  {q{x,-.x,)ln(X,/X2) 

JC,X|(X, -f  X|) 

— 5^ — 

(^2(x,x)  =  (/)  +q){X2-X2){0,  +  e2) 

+  l9ix2-X2)]n(X,IX2) 

XiXtiXi-Xi) 


-Ip 


1 


^  47t(l— v) 

0i  =  arctan 


^2 

,  9=;2(1-2v), 


>  (10) 


Vi+^J’  2 


<  01  ct 


<02  <2 


>■  (11) 


02  =  arctan 

X,  =  (X,  -X,y  +  (X2  -  X2)^ 
^2  =  (x,  +  X,  -F  (x2  - 


To  solve  equation  (8)  for  Fj,  we  write  <^j(x,  x)  a 
GM3(x,x).The  stress  is  then 

^i3(x.X)  =  ^:  =  (12 

All  other  stress  components  are  zero.  For  a  uni 
concentrated  force  Fj  at  x,  equation  (8)  gives 


V^d'3  +  ^(x-x)  =  0  (13) 

where  S(x,  x)  is  the  Dirac  Delta  function.  With  the 
boundary  conditions  of 

^<l>^  . 

-r —  =  0  at  X|  =  0  and 

503  507  ^  - 

—  =  0  and  — -->0  as  +  (14) 

we  have 

03(x,  x)  =  (In  X,  +  In  X2)  (15) 

where  ,  X2  are  given  in  equation  (1 1).  Hence  the 
stress  fields  caused  by  the  equivalent  forces  due  to  slip 
strain  can  readily  be  calculated. 


4.  EFFECT  OF  THE  INHIBITION  OF  CROSS-SLIP 
OR  MULTIPLE  SLIP  IN  INTERMETALLIC 
COMPOUNDS  ON  HIGH-CYCLE  FATIGUE 

The  above  gives  a  method  of  calculation  of  the 
growth  of  an  extrusion  or  an  intrusion  due  to  a  single 
slip  system.  To  understand  the  effect  of  the  inhibition 
of  cross  or  multiple  slip  of  ordered  intermetallic 
alloys  on  fatigue  crack  initiation,  we  need  to  study 
the  effect  of  secondary  slip  on  the  extent  of  extrusion. 


The  buildup  of  this  local  slip  strain  e^p  in  P  and  Q 
is  caused  by  the  positive  and  negative  initial  shear 
stress  x\p  which,  in  turn,  is  caused  by  e\^  in  R. 
If  R  were  cut  out,  this  free  length  of  R  would  be 
longer  than  the  slot  cut,  by  an  amount  referred  to  as 
the  “static  extrusion”  [1 1].  This  causes  an  initial 
compression  in  R,  As  the  extrusion  grows  under 
cyclic  loading,  the  slice  R  increases  in  length.  This 
elongation  causes  the  compression  to  decrease.  A 
question  has  been  raised  as  to  whether  the  extrusion 
growth  will  cease  after  the  extrusion  has  reached  the 
static  extrusion.  The  residual  tensile  stress  caused 
by  elongation  in  R  due  to  extrusion  can  cause  changes 
of  resolved  shear  stresses  in  all  12  slip  systems.  The 
resolved  shear  stress  in  one  slip  system  may  reach  the 
critical  shear  stress  and  slide.  The  plastic  strain 
caused  by  slip  in  this  system  has  a  tensor  component 
just  like  e\ai  in  causing  the  positive  and  negative 
respectively,  in  P  and  Q.  This  secondary  slip  has 
been  shown  [16, 17]  to  increase  greatly  the  extent  of 
extrusion  and  intrusion.  This  method  is  here  used  to 
calculate  the  fatigue  damage  for  the  case  with  second 
slip  and  that  without  second  slip.  These  two  cases 
correspond  to  disordered  and  ordered  superalloys 
respectively. 


5.  NUMERICAL  CALCULATION 


The  mechanical  properties  used  in  the  calculation 
were  taken  to  correspond  to  NigAI.  The  critical  re¬ 
solved  shear  stress  t*'  was  taken  to  be  200  MPa,  Pope 
and  Ezz  [18].  The  shear  modulus  G  was  calculated 
from  the  elastic  compliance  5,2,2  [18].  This  5,2,2  is 
0.848  X  10“^  (1/GPa),  gives 


^*^1212 


10^  MPa 

4  X  0.848  X  10“^ 


=  29,400  MPa. 


The  Possion’s  ratio  was  taken  to  be  0.3.  Calculations 
have  been  made  for  a  single  fatigue  band  and  then 
multiple  fatigue  bands.  Since  the  local  plastic  strain 
in  the  band  is  much  larger  than  the  average  plastic 
strain  in  the  crystal,  the  strain-hardening  rate  in 
the  band  is  much  less  than  that  of  the  crystal.  The 
strain-hardening  in  the  slip  band  is  here  neglected. 


(a)  A  single  fatigue  band 

The  dimensions  of  this  single  fatigue  band  in  the 
most  favorably  oriented  crystal  at  the  free  surface 
of  an  intermetallic  compound  polycrystal  are  shown 
in  Fig.  1.  The  initial  resolved  shear  stress  t’  was 
taken  to  be  0.46  MPa.  =  200.0  MPa.  This  gives 
an  excessive  shear  stress  t^,  defined  as  t®  -I-  t’  ~  of 
0.46  MPa.  The  three  slip  directions  of  each  of  the 
four  slip  planes  of  a  f.c.c.  crystal  are  shown  in  Fig.  2. 
The  most  favorably  oriented  crystal  of  a  polycrystal 
loaded  under  alternate  tension  and  compression 
along  the  Xj-axis  (Fig.  1)  has  a  slip  plane  and  a  slip 
direction  making  an  angle  of  45^  with  the  direction  of 
loading.  The  slip  system  is  referred  to  as  the  primary 
slip  system.  Let  a,  in  Fig.  2,  correspond  to  this  system. 
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Fig.  2.  Crystallographic  direction  of  a  f.c.c.  crystal. 


During  fatigue  loading,  the  build-up  of  large  local 
plastic  shear  strain  in  the  primary  slip  system,  positive 
in  P  and  negative  in  Q,  tends  to  start  an  extrusion 
or  an  intrusion.  Consequently  an  appreciable  direct 
stress  r„  will  be  built  up  in  R.  The  Schmid  factors 
of  all  the  12  slip  systems  for  the  stress  r.y  and  those 
for  have  been  calculated  [16]  as  shown  in  Table  1 . 
It  was  found  that  there  are  four  slip  systems  C| ,  c^, 
</i  and  equally  favorably  under  t„.  Of  these 


four,  C3  is  most  favorably  under  T32,  where  T22  is 
the  alternate  loading.  Hence  c,  is  considered  to  be  the 
active  second  slip  system  in  R. 

For  numerical  calculations,  the  slices  P,  Q,  R  are 
divided  into  a  number  of  parallelogram  grids  along 
the  length  in  the  a-direction.  Each  grid  is  taken  to 
have  constant  plastic  strain.  The  average  stress  in  the 
grid  is  taken  as  the  corresponding  grid  stress.  From 
the  plane  strain  solution  of  a  semi-infinite  medium  the 
stress  field  caused  by  a  uniform  plastic  strain  in 
the  «th  grid  was  calculated.  The  average  relieved 
stress  over  the  mth  grid  denoted  by  is 

written  as 

=  ~C{m,  xP;  n,  a.p)el^^. 

As  a  cyclic  loading  proceeds,  slips  occur  in  P  and 
Q  and  cause  a  residual  stress  field.  The  total  resolved 
shear  stress  In  the  mth  grid  of  xp  slip  system  is  the 
sum  of  the  initial,  residual  and  applied  stresses 

-  Y.  C(m,  xp-,n,  (16) 

It 

where  n  is  summed  over  all  the  sliding  grids  with 
plastic  strain 

In  the  grids  where  t(jc)=+t',  an  incremental 
applied  shear  stress  At“  ,  with  At„  =  At  =  0,  yields 

=  -I  C{m  xP;n,  xp)  (17) 


Table  1.  Resolved  shear  stresses  in  different  slip  systems  caused  by  cyclic  loading  tjj  and  t„ 
_ _ (g  denotes  02  direction)  ” 


Slip 

system 

Norma!  to  slip  plane 
^(m) 

Slip  direction 
n{m) 

Schmid  factor 

0\ 

{4-4  4) 

0 

(##«) 

0 

0 

h, 

{44-4) 

(41  42  42  42  42\ 

V  4  ^  3  ’  4  3  ’  ~) 

0 

*2 

_2n/2  ^ 

^3^,^,3^,+42  43'^ 

0 

A3 

{44‘) 

0 

‘•j 

(^.s/2  43  42 
\3  *  6'  3  6 '  3) 

{444) 

0.408 

^■2 

(41.42  42  42  42\ 

U  6  ’  3  6  •  3) 

0 

(4^  ,42  42  42  \/2\ 

1  3  6  ’  3  6  ’  ~—) 

(-41+42  42  42  42\ 

\  4  3  •"  4  •  3  '~) 

-0.408 

^1 

(42  42  s/2_^42  4a 

\3  6  ■  3  +  6  '  6  j 

(-41  42  42.42  n/3\ 

\  4  3-4 

-0.408 

d2 

( 42.  _  42  42  42  4^ 

1  3  6  ’  3  +  6  ’Tj 

{444) 

0.408 

dy 

(42  42  42 ,42  42\ 

V  3  6-3  — j 

{444) 

0 
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There  are  as  many  unknown  as  the  number  of 
the  above  equations.  The  plastic  strain  increments 
in  the  sliding  grids  for  an  incremental  loading 
can  be  readily  calculated.  Similarly,  the  stress 
components  and  Xpp  in  R  can  be  calculated  [14] 
and  written  as 


— =  L  Cim,  aa;  n,  ' 

n 

-■ =  E  <^(w,  n, 


(18) 


As  discussed  earlier,  is  the  most  likely  second  slip 
system  to  become  active.  The  resolved  shear  stress  in 
Cj  is  denoted  by  The  initial  resolved  shear  stress 
T'f,  varies  with  the  initial  stress  components  t\3, 
t*33.  For  our  numerical  calculations,  is  assumed  to 
be  zero.  The  incremental  resolved  shear  stress  in  this 
second  slip  system 


=  0.31  8At22  -  0.408At^  .  (19) 

When  increases  to  t*'  or  decreases  to  —  this  slip 
system  slides  and  causes  Now  we  have  slip  in  the 
(xp  slip  system  in  P  and  Q  and  in  slip  system  in  R. 

Let  the  two  slices  P  and  Q  be  divided  into  2N  grids 
and  R  into  M  grids.  The  resolved  shear  stress  in  the 
up  and  the  slip  system  of  the  /th  grid  is  written  as 


=  -C{i,  ap;j\  -  ca,  ap;  k, 

(^PMpj  -  CO,  k, 


(20) 


The  repetition  of  subscript  J  denotes  summation  from 
1  to  2N  and  that  of  subscript  k  denotes  summation 
from  1  to  M.  C(/,  aP)  is  the  residual  stress  tJ,, 
in  the  /th  grid  due  to  a  unit  plastic  strain  in  the 
yth  grid.  The  incremental  plastic  strain  in  all  grids  are 
calculated  from  equation  (17). 

With  the  method  described,  the  extrusion  growth 
of  a  single  fatigue  band  with  cycles  of  loading  was 
calculated  for  the  case  with  no  slip  in  the  second  slip 
system,  corresponding  to  an  ordered  superalloy  and 
that  with  slip  in  the  second  slip  system  corresponding 
to  a  disordered  compound.  These  results  are  shown  in 
Fig.  3.  It  is  seen  that  ordered  intermetallic  compound 
has  a  much  lower  rate  of  extrusion  on  intrusion 
growth. 


Fig.  3.  Plastic  shear  strain  at  free  surface  vs  number  of 
cycles,  single  fatigue  band. 


^22 


Grain 

boundary 


Fig.  4.  Spacing  of  fatigue  bands  and  their  initial  stress. 


(b)  Multiple  fatigue  bands 

A  number  of  slip  bands  have  been  observed  in 
a  most  favorably  oriented  crystal  in  a  polycrystal 
subject  to  cyclic  loading  [5],  The  crystals  surrounding 
this  most  favorably  oriented  crystal  are  less  favorably 
oriented  and  hence  have  lower  resolved  shear  stresses 
caused  by  the  applied  load.  In  the  early  stages  of  a 
high-cycle  fatigue,  the  residual  resolved  shear  stress 
in  the  neighboring  crystals  is  small  compared  to  the 
decrease  of  the  applied  resolved  shear  stress  due  to  the 
difference  in  orientation.  Hence,  in  the  early  stages  of 
high-cycle  fatigue,  slip  in  the  most  favorably  oriented 
crystal  can  develope  to  some  extent  without  signifi¬ 
cant  slip  in  the  neighboring  crystals.  For  numerical 
study,  we  consider  the  most  favorable  crystal  to  have 
eleven  slip  bands  spaced  at  3  microns  apart  as  shown 
in  Fig.  4.  Each  band  is  assumed  to  consist  of  three 
thin  slices  P,  R  as  shown  in  Fig.  1 .  P  and  Q 

in  each  band  are  assumed  to  have  equal  and  opposite 
uniform  initial  resolved  shear  stresses  of  the  amounts 
as  indicated  in  Fig.  4.  The  initial  excessive  shear  stress 
T^isT‘4-T®  —  T^.  During  loading  this  excessive  shear 
stress  is  to  be  relieved  with  slip  in  the  different  slip 
bands.  Slip  distributions  in  these  bands  are  readily 
calculated.  Slip  in  one  band  not  only  relieves  shear 
stress  in  this  band  but  also  in  the  other  bands.  This 
gives  the  interaction  of  slip  bands. 

Using  the  method  described  above,  each  thin  slice 
P,  Q  and  R  of  each  band  is  divided  along  the  a 
direction  into  three  parallelogram  grids,  and  the  slip 
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Number  of  cycles 

Fig.  5.  Plastic  strain  at  free  surface  vs  number  of  cycles  at  different  applied  excessive  shear  stresses. 


distributions  in  the-  different  slices  of  the  slip  bands  shear  stresses  for  both  ordered  and  disordered  inter- 

were  calculated.  The  center  fatigue  band  has  the  metallic  compounds  or  alloys  as  shown  in  Fig.  5.  For 

maximum  initial  shear  stress  and  has  been  found  to  each  value  of  this  excessive  shear  stress,  the  number 

give  the  maximum  plastic  shear  strain  at  the  free  of  active  slip  bands  at  different  cycles  of  loading  are 

surface.  This  shear  strain  is  considered  to  be  a  shown  in  Fig.  6.  From  the  curves  in  Fig.  5  we  pick 

measure  of  the  fatigue  damage.  The  variation  of  this  the  number  of  cycles  to  yield  the  given  value  of  the 

surface  strain  with  the  number  of  cycles  of  loading  plastic  resolved  shear  strain,  e"p  at  the  free  surface, 

are  plotted  at  different  values  of  excessive  resolved  thus  we  pbtain  the  excessive  shear  stress  vs  cycles 


N 


400  800  1200  1600  2000 

Number  of  cycles 

Fig.  6.  Number  of  active  fatigue  bands  at  free  surface  vs  number  of  cycles. 
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Fig.  7.  Maximum  initial  excessive  shear  stress  vs  number 
of  cycles  curves  at  different  levels  of  fatigue  damage. 

of  loading  to  yield  the  given  surface  plastic  strain, 
which  is  a  measure  of  the  extent  of  extrusion  or 
intrusion  and  is  taken  as  the  extent  of  crack  initiation. 
These  vs  Af  curves  for  both  the  ordered  and 
disordered  alloys  are  shown  in  Fig.  7.  It  is  seen 
that  the  shap>e  of  these  curves  for  different  levels 
of  fatigue  damage  is  somewhat  similar  to  a  typical 
S-N  curve. 

At  the  initial  stage  of  fatigue  loading,  the  slip  in  a 
single  band  is  directly  proportional  to  the  excessive 
shear  stress.  The  rate  of  incremental  slip  in  P  and 
Q  per  cycle  depends  on  the  positive  and  negative 
resolved  shear  stresses  (t‘  +  As  loading  proceeds, 
extrusion  grows  and  this  shear  stress  (t*  +  tO 
decreases  and  reduces  the  growth  rate  of  extrusion. 
This  has  the  effect  of  decreasing  the  slope  of  the  e"p 
at  the  free  surface  with  N.  It  has  been  shown  that  this 
slope  approaches  zero  when  the  extrusion  approaches 
the  “static  extrusion”  [17].  When  more  than  one  slip 
band  becomes  active,  slip  in  one  band  relieves  the 
shear  stress  not  only  in  this  band  but  also  in  other 
bands.  This  decreases  the  amount  of  slip  in  each 
band.  This  causes  the  curves  e"p  vs  N  to  bend  further 
down.  This  has  the  effect  of  causing  the  vs  N  curve 
to  give  a  specified  value  of  surface  strain  to  become 
more  flat. 

Stress  amplitude  vs  cycles  to  fatigue  failure  from 
tests  by  Boettner  et  al.  [3]  for  ordered  and  quenched 
Ni3Mn  at  25°C  is  shown  in  Fig.  8.  These  test  curves 
of  Ni3  Mn  instead  of  these  of  Nij  Al  are  here  used  to 
compare  the  calculated  results,  since  the  correspond¬ 
ing  test  data  of  Ni3AI  is  not  found  in  the  literature. 
The  trend  of  these  tests  clearly  agree  to  the  calculated 
results  as  shown  in  Fig.  7.  This  explains  why  an 
ordered  intermetallic  compound  or  alloy  has  longer 


fatigue  life  than  the  disordered  ones.  Since  the  yield 
stress  of  ordered  intermetallic  alloys  is  lower  than  the 
disordered,  the  longer  fatigue  life  is  due  to  the  longer 
life  of  crack  initiation.  This  seems  to  explain  well  this 
observed  longer  high-cycle  fatigue  life  of  the  ordered 
intermetallic  compounds. 
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Abstract 


Tests  have  shown  that  ordered  intermetallic  compounds  favor  planar  slip  and  tend  to  have  no  cross-slip  or  multislip.  The 
yield  stress  of  ordered  compounds  is  less  than  that  of  disordered  compounds.  However,  the  stress-controlled,  high  cycle 
fatigue  life  of  ordered  compounds  is  appreciably  more  than  that  of  disordered  compounds.  This  is  considered  to  be 
caused  by  a  low  rate  of  crack  initiation.  In  contrast,  strain-controlled  tests  show  that  ordering  shortens  life  by  crack  initia¬ 
tion.  In  this  paper,  we  attempt  to  explain  these  observed  behaviors  by  the  micromechanic  theory  of  high  cycle  fatigue 
crack  initiation  developed  by  Lin  and  coworkers.  The  effect  of  ordering  on  crack  initiation  under  cyclic  combined  axial 
and  torsional  loadings  is  also  shown. 


1.  Introduction 

Ordered  intermetallic  compounds  are  distinct  from 
ordinary  alloys  according  to  their  lattice  structure  and 
bonding.  If  the  atoms  of  the  two  elements  of  the  com¬ 
pounds  segregate  fairly  completely  to  designated 
atomic  sites,  then  it  is  called  an  ordered  solution  or 
superlattice. 

At  low  and  intermediate  temperatures,  nickel-based 
intermetallic  compounds  with  LI 2  structures  deform 
by  (1 10){1 1 1}  slip.  If  the  slip  occurs  by  the  motion  of 
KllO){llll  dislocations,  they  will  leave  a  plane  of 
antiphase  boundary  in  their  wake.  The  energy  required 
to  create  such  a  fault  is  high,  which  makes  this  slip 
difficult.  In  contrast,  slip  by  the  motion  of  (110>{111} 
dislocations  does  not  create  such  a  fault,  so  occurs 
frequently.  However,  because  the  line  energy  of  a 
dislocation  is  proportional  to  the  square  of  the  distance 
which  the  atom  moves,  this  superdislocation 
(1 10){1 1 1 }  needs  high  energy  to  move.  Nevertheless,  if 
the  superdislocation  dissociates  into  two  superpartial 
dislocations,  this  dislocation  energy  is  reduced.  There¬ 
fore,  the  energy  of  a  superdislocation  depends  on  the 
length  of  the  Burgers  vector  and  the  antiphase  bound¬ 
ary  created  by  slip.  Dislocations  in  a  superlattice  tend 
to  move  in  a  group,  connected  by  a  strip  of  antiphase 
boundary.  The  leading  dislocation  creates  an  antiphase 
boundary  and  the  trailing  dislocation  eliminates  it.  It  is 
unlikely  that  cross-slip  of  the  associated  pair  can  occur, 
since  this  requires  the  trailing  dislocation  to  follow 
exactly  in  the  wake  of  the  leading  dislocation.  This 


causes  cross-slip  or  multislip  to  be  difficult  in  ordered 
intermetallic  compounds. 

Nickel-base  intermetallic  compounds  with  a  long- 
range  ordered  structure  have  been  found  to  exhibit  an 
anomalous  increase  in  flow  stress  with  temperature.  At 
a  low  temperature,  these  compounds  deform  by 
(1 10){  1 1 1 }  slip;  however,  at  elevated  temperatures,  slip 
occurred  on  the  {010}  plane.  Various  theories  have 
been  advanced  to  explain  the  observed  temperature 
dependence  of  the  yield  stress.  Copley  and  Kear  [1] 
proposed  that,  at  low  temperatures,  the  superpartial 
dislocations  are  dissocaited  to  Shockley  partid  disloca¬ 
tions,  such  that  flow  occurs  by  the  easy  i(112){lll} 
shear  of  the  crystal.  At  elevated  temperatures,  the 
superpartials  are  much  less  widely  dissociated,  leading 
to  the  difficult  1(110){111}  shear  of  the  crystal  and, 
hence,  giving  a  high  yield  strength. 


2.  Experimental  results 

2.1.  Stress-controlled  fatigue 
Tests  have  shown  that  long-range  order  substantially 
increases  the  high  cycle  fatigue  lives  of  FeCo-V,  Ni3Mn 
(Fig.  1)  and  other  intermetallic  compounds,  despite  a 
decrease  in  yield  stress  [2].  The  ratio  of  the  fatigue  limit 
to  the  yield  strength  exceeds  0.5  for  most  of  the 
ordered  alloys;  in  contrast,  under  the  disordered  con- 
diion,  the  ratio  decreases  significantly.  Since  crack 
propagation  is  much  more  rapid  in  the  ordered  con- 
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dition,  it  has  been  concluded  by  Stoloff  et  al.  [3]  that 
ordering  substantially  delays  crack  initiation. 

2.2.  Strain-controlled  fatigue 

Hsiung  and  Stoloff  [4]  have  conducted  experiments 
on  fatigue  crack  initiation  in  two  intermetallic  com¬ 
pounds— NijAl  single  crystals  and  FeCo-V  polycrys¬ 
tal— under  strain-controlled  cycling.  Strain-controlled 
tension-compression  cyclic  deformation  tests,  begin¬ 
ning  in  tension,  were  performed  on  a  closed  lop 
machine.  The  tests  showed  that  order  substantially 
decreases  the  strain-controlled  fatigue  resistance  in  this 
alloy  (Fig.  2).  However,  there  is  little  difference  in  the 
development  of  dislocation  substructures  and  crack 
initiation  sites  [4].  Bundles  of  unit  dislocations  oriented 
crystallographically  were  found  in  the  disordered 
material.  Dislocation  pairs  were  seen  in  the  ordered 
alloy  samples.  Slip  along  (110)  occurs  in  both  ordered 
and  disordered  FfeCo-V,  and  cracks  initiated  almost 
exclusively  at  slip  bands,  intrusions  and  extrusions. 
Tests  carried  out  in  high  vacuum  also  have  shown  that 
slip  bands  are  preferred  sites  for  crack  initiation  and 
propagation  in  both  the  ordered  and  disordered  condi¬ 
tions.  The  height  of  extrusion  and  the  depth  of  intru¬ 
sion  are  taken  as  measures  of  the  fatigue  crack 
initiation. 

In  short,  the  test  results  show  that  long-range 
ordered  substantially  increases  the  fatigue  lives  of 
intermetallic  compounds  in  stress-controlled  cycling 
but  decreases  the  strain-controlled  fatigue  resistance 
(Figs.  1  and  2).  Stoloff  et  al.  [3]  indicated  that  the  least 
understood  aspects  of  the  cyclic  behavior  of  inter- 


Fig.  1.  Stress  amplitude  or  cycles  to  failure  for  Ni,Mn  fatigued  at 
25  "C  [2]. 


metallic  compounds  are  the  processes  leading  to  crack 
initiation.  This  paper  is  an  attempt  to  explain  the 
observed  high  cycle  fatigue  crack  initiation  character¬ 
istics  of  these  intermetallic  compounds. 


3.  Micromechanic  quantitative  theory  of  fatigue  crack 
initiation  [5, 6] 

Extrusion  and  intrusion  were  first  discovered  by 
Forsyth  and  Stubbington  in  1955  [7].  For  an  extrusion 
to  start,  positive  shear  has  to  occur  in  a  thin  slice  (P) 
and  negative  shear  in  a  closely  located  slice  (Q),  as 
shown  in  Fig.  3.  Initial  defects  exist  in  all  metals  and 
cause  an  initial  stress  field  r'.  An  initial  shear  stress 
field  favorable  to  positive  shear  in  P  and  negative  shear 
in  Q  is  clearly  one  having  a  positive  shear  stress  in  P 
and  a  ngative  shear  stress  in  Q.  Such  an  initial  stress 
field  can  be  provided  by  an  initial  tensile  strain  in 
R  (The  repetition  of  Greek  letters  does  not  denote 
summation.) 

Let  us  consider  a  perfect  crystal.  If  we  cut  a  slit 
through  a  crystal  and  force  a  sheet  of  metal  one  atom 
thick  into  the  slit,  a  pair  of  parallel  edge  dislocations  (A 
and  B)  of  opposite  signs,  forming  an  interstitial  dipole, 
is  produced,  as  shown  in  Fig.  4.  If  we  cut  a  rectangular 
block  along  the  dotted  line,  the  free  length  of  this  block 
will  be  one  atomic  spacing  more  than  the  correspond¬ 
ing  length  of  the  hole.  If  there  are  n  such  dipoles  in  a 
length  of  atomic  spacing  N,  this  will  give  an  initial 
strain  of  n/N.  Hence,  this  initial  strain  can  be 
caused  by  an  array  of  dislocation  dipoles.  This  array  of 
dipoles  was  suggested  by  Lin  and  Ito  [8]  in  1969  as  a 
possible  way  of  providing  the  initial  strain  to  cause  the 
favorable  initial  stress  field. 

It  follows  that  this  positive  can  be  caused  by  a 
row  of  interstitial  dislocation  dipoles  and  a  negative  ei,„ 
by  vacancy  dipoles  [8].  It  was  shown  by  Mughrabi  et  al. 
[9]  that  the  experimentally  observed  ladder  structure  in 
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Fig.  3.  Schematic  diagram  of  most  favorably  oriented  crystal  at 
free  surface. 
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Fig.  4.  Schematic  diagram  of  a  disloation  interstitial  dipole. 


the  persistent  slip  bands  (PSBs)  indicates  dislocation 
dipoles.  These  cause  such  an  initial  resolved  shear 
stresses  at  the  interface  between  the  PSB  and  the 
matrix. 

If  a  metal  specimen  is  loaded  in  cyclic  tension  and 
compression  of  low  amplitude,  the  plastic  deformation 
essentially  occurs  only  in  the  most  favorably  oriented 
crystal.  This  crystal  is  considered  to  be  at  a  free 
surface,  with  a)3  as  the  most  favorably  oriented  slip 
system,  as  shown  in  Fig.  3.  During  loading,  when  the 
resolved  shear  stress  in  some  region  reaches  the  critical 
shear  stress  r',  slip  occurs.  After  loading,  this  slip 
remains  and  causes  a  residual  resolved  shear  stress  t\ 


Then,  we  have  the  total  shear  stress  equal  to  the  sum  of 
these  stresses,  such  that 

r=r'  +  T^+r'  (1) 

where  r’  is  the  applied  shear  stress. 

3. 1.  Monotonic  building  of  extrusion  and  intrusion 
15,6/ 

Tests  show  that  the  extrusions  and  intrusions 
increase  monotonically  with  cyclic  loading.  For  such  an 
extrusion  to  build-up  monotonically,  slip  in  P  (Fig.  3) 
during  forward  loading  has  to  be  irreversible  during  the 
reversed  loading.  This  clearly  needs  a  gating  mecha¬ 
nism  to  prevent  slip  from  reversing.  It  is  interesting  to 
see  that  the  microstress  field  produced  by  slip  supplies 
this  gating  mechanism. 

Let  us  consider  an  initial  stress  field  r'  to  be  positive 
in  P  and  negative  in  Q,  as  explained  previously.  A  posi¬ 
tive  loading  T22  on  the  polycrystal  (Fig.  3)  causes  a 
positive  r“  in  the  whole  crystal,  r'  +  r®  reaches  the 
critical  shear  stress  r®  first  in  P  and  so  this  slides  first. 
The  plastic  strain  resulting  from  this  slip  is  taken  to  be 
constant  along  the  axis  (Fig.  5).  This  gives 


The  equilibrium  condition  r,yy=0  in  rectangular 
coordinates  x,,  X2  and  X3  in  terms  of  coordinates  a,  P 
and  X3  gives 

3a  3^ 

The  first  term  here  is  finite.  Since  the  sum  of  the  two 
terms  vanishes,  the  second  term  also  must  be  finite. 
Hence,  the  change  in  along  )3,  across  the  small 
distance  from  P  to  Q,  owing  to  slip,  is  very  small  [5, 6]. 

Therefore,  the  slip  in  P  relieves  not  only  the  positive 
shear  stress  in  P  but  also  that  in  its  neighboring  region 
including  Q.  This  slip  increases  the  negative  resolved 
shear  stress  in  Q,  causing  Q  to  slide  more  readily 
during  the  reverse  loading  [10].  Similarly,  the  negative 
slip  in  Q  relieves  the  negative  shear  stress  not  only  in  Q 
but  also  in  P,  thus  causing  P  to  slide  more  readily 
during  the  second  forward  loading.  This  process  is 
repeated  for  every  cycle,  thus  providing  a  natural 
gating  mechanism  to  cause  the  monotonic  build-up  of 
the  local  slip  strain  in  P  and  Q,  as  observed  in 
experiments,  pushing  R  out  of  the  free  surface  to  start 
an  extrusion.  Interchanging  the  signs  of  the  initial 
stresses  initiates  an  intrusion  instead  of  an  extrusion. 
This  gating  mechanism  is  briefly  reviewed  here  for  the 
readers’  convenience.  This  gating  mechanism  was 
introduced  by  Lin  and  Ito  [10]  in  1969,  was  expressed 
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Fig.  5.  Crystallographic  directions  of  an  f.c.c.  crystal. 


by  dislocation  arrays  by  Tanaka  and  Mura  [1 1  ]  in  198 1 
and  is  also  discussed  in  Suresh's  excellent  book  Fatigue 
of  Materials  [12].  This  mechanism  depends  on  the 
closeness  of  slid  slices  P  and  Q.  This  explains  why  the 
extrusions  are  very  thin,  as  observed  in  fatigue  speci¬ 
mens. 

The  governing  condition  for  the  initiation  or  con¬ 
tinuation  of  slip  is  having  the  resolved  shear  stress  t 
equal  to  the  critical  shear  stress  r^.  The  residual  stress 

is  a  component  of  the  resolved  shear  stress  t.  To 
calculate  the  analogy  of  plastic  strain  and  applied 
force  developed  by  Lin  [13]  is  used.  Plastic  strain  is 
generally  caused  by  the  displacement  of  dislocations. 
As  indicated  by  Hirth  and  Lothe  in  their  excellent 
book  Theory  of  Dislocations  [14],  most  dislocation 
literature  deals  with  dislocations  in  the  isotropic 
approximation.  They  have  given  two  good  reasons  for 
these  approximations.  “First,  the  mathematics  becomes 
much  more  complicated  when  anisotropy  is  consid¬ 
ered.  Second,  in  many  cases  the  errors  involved  in  the 
isotropic  approximation  are  about  20  to  30%,  and 
these  are  submerged  by  other  approximations  in  dis¬ 
location  theory  or  by  errors  in  experimental  observa¬ 
tion.” 

Hence,  this  isotropic  approximation  is  not  likely  to 
change  the  general  results,  as  shown  in  Figs.  6  and  7. 
To  simplify  the  numerical  calculation,  this  anisotropy 
of  the  intermetallic  is  neglected.  Then,  the  residual 
stress  resulting  from  a  given  plastic  strain  distribution 
e]i  is  given  by 

n  =  +  2  ~  U,/,,  -IGe^  (2) 

where  is  the  strain  caused  by  d^j  is  the 
Kronecker  delta,  and  A  and  G  are  Lame’s  constants. 
The  length  of  slices  P,  Q  and  R  along  the  direction 


Humber  of  Cycles 

Fig.  6.  Maximum  initial  excessive  shear  stress  vs,  number  of 
cycles,  Le,  t-N  curves,  at  different  levels  of  fatigue  damage. 


Fig.  7.  Variation  of  cycle  stress  with  cycles  of  loading  in  strain- 
controlled  fatigue. 

are  much  larger  than  the  thickness  and  their  inclined 
lengths  shown  in  Fig.  3. 

Strains  in  the  major  central  portion  of  the  slices  are 
taken  to  be  of  generalized  plane  strain.  The  displace- 
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ment  along  the  x.  axis  is  denoted  as 


4.  Calculated  results 


=  as)  ;=  1,2,3  (3) 

The  stress-strain  relationship  is  then  [5, 6] 

2G  j-pr^  +  2 

where  6=  u^^  +U2_2  and  v  is  Poisson's  ratio. 

Substituting  eqn.  (4)  into  the  condition  of  equilib¬ 
rium  yields 


v-«,+ 


1  de 

l-2vdXj 


v-«, 


0 


/=1,2 


(5) 

(6) 


where 


(7) 


Equations  (5)  and  (6)  are  not  coupled  and  can  be 
solved  separately.  The  details  of  these  solutions  are 
given  in  refs.  5  and  6. 


3.2.  Effect  of  inhibition  of  cross-slip  or  multiple  slip  on 
fatigue  crack 

As  indicated  previously,  cross-slip  or  multiple  slip  is  . 
inhibited  in  ordered  intermetallic  compounds.  The 
build-up  of  local  plastic  strain  in  P  and  Q  depends 
on  the  initial  stress  r',  which  is  caused  by  e'^a  in  R.  This 
e'oa  is  due  to  the  extra  initial  length  of  R  over  that 
provided  by  the  slot  between  P  and  Q.  An  extrusion 
grows  at  the  expense  of  this  extra  length.  When  the 
extrusion  equals  this  extra  length— which  is  then  called 
a  static  extrusion  [9] — the  growth  may  stop.  However, 
the  growth  of  an  extrusion  causes  a  residual  tensile 
.stress  in  R.  This  can  cause  changes  in  the  resolved 
shear  stress  in  all  12  slip  systems  of  the  f.c.c.  cry.stal  [5, 
6].  The  resolved  shear  stress  in  one  slip  system  may 
reach  the  critical  shear  stress  and  slide.  The  plastic 
strain  caused  by  this  slip  has  a  tensor  component 
e'aa,  which  is  similar  to  e'^a  in  causing  the  po.sitive  and 
negative  in  P  and  Q  respectively.  This  e'j,  has  been 
shown  to  increase  greatly  the  extent  of  an  extrusion 
[5,6]. 

The  above  gives  the  analytical  method  to  calculate 
the  growth  of  an  extrusion  and  intrusion  for  the  cases 
with  secondary  slip  (corresponding  to  disordered)  and 
without  secondary  .slip  (corresponding  to  ordered). 
These  calculations  were  carried  out  for  both  stress- 
controlled  and  strain-controlled  loadings. 


The  mechanical  properties  used  in  the  present 
calculations  were  taken  to  corre.spond  to  Ni^AI  [15]. 
The  critical  resolved  shear  stress  t'  was  taken  to  be 
200  MPa  [16].  The  .shear  modulus  G  was  obtained 
from  the  mea.sured  elastic  compliance  5,212  17, 

giving  G=29  400  MPa.  Poisson’s  ratio  v  was  taken  to 
be  0.3.  The  strain-hardening  rate  in  the  band  is  much 
less  than  that  in  the  crystal.  Hence,  the  .strain  hardening 
in  the  slip  band  is  neglected. 

The  dimensions  of  the  fatigue  band  in  the  most 
favorably  oriented  crystal  at  the  free  surface  of  an 
intermetallic  compound  polycrystal  are  shown  in  Fig. 
3.  The  initial  resolved  shear  stress  r'  in  P  and  Q  was 
taken  to  be  ±0.46  MPa.  The  applied  shear  stress 
r“  =  200  MPa.  This  gives  an  excessive  shear  stress 
defined  as  r®  +  r’  —  of  0.46  MPa.  The  slip  directions 

and  slip  planes  of  an  f.c.c.  crystal  are  indicated  in  Fig.  5. 
The  most  favorably  oriented  crystal  of  the  polycrystal 
has  a  slip  system  with  a  slip  direction  and  slip  plane 
making  an  angle  of  45°  with  the  loading  axis  (Fig.  3). 
This  slip  system  is  identified  as  02  in  Fig.  5  and  is 
referred  to  as  the  primary  slip  system.  During  fatigue 
loading,  an  extrusion  grows,  causing  a  tensile  stress  in 
R.  The  Schmid  factors  of  all  12  slip  systems  for  the 
loadings  r22  and  x^a  have  been  calculated  [5,  6].  It  was 
found  that  four  slip  systems,  i.e.  c,,  Cj,  tf,  and  d2,  are 
equally  favorable  under  x^^.  Of  these  four,  C3  is  the 
most  favorable  under  T22.  Hence,  Cj  is  taken  to  be  the 
active  second  slip  system  in  R. 

The  slices  P,  Q  and  R  are  divided  into  a  number  of 
parallelogram  grids  along  the  length  in  the  a  direction. 
A  constant  plastic  strain  is  assumed  in  each  grid.  The 
average  stress  in  the  grid  is  taken  as  the  corresponding 
grid  stress.  From  the  plane  strain  solution  of  a  semi¬ 
infinite  medium,  the  stress  field  caused  by  the  uniform 
plastic  strain  in  the  nth  grid  was  calculated.  The 
average  relieved  stress  over  the  mth  grid  is  expressed 
as 

'^np„=~C{m,a^;n,aP)e^f,^  (8) 

Here,  the  coefficient  C(m,  aff,  n,  afi),  for  the  stress 
Xap^  relieved  in  the  mth  grid  owing  to  a  unit  shear 
strain  equals  C  (n,  ap-,  m,  aP)  for  the  stress 
relieved  in  the  nth  grid  owing  to  a  unit  shear  strain 
This  reciprocal  relationship  facilitates  numerical 
calculation. 

In  the  grids  where  x(x)  =  an  incremental  applied 
shear  stress  is  then 

=  -S.C{m,  aP\  n,  aP)Le"„p^  (9) 

This  gives  a  set  of  equations  to  solve  the  unknown 

values.  The  initial  resolved  shear  stress  in  Cj, 
denoted  by  r'j,,  is  assumed  to  be  zero.  For  each 
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unknown  there  is  one  of  the  above  equations. 

Hence,  the  plastic  strain  increments  in  the  sliding  grids 
*  for  an  incremental  applied  load  are  readily  calculated. 
Indeed,  the  incremental  shear  stress  in  the  second  slip 
system  Cj  is  expressed  as 

Ar^,  =  0.318AT22-0.408AT„„  (10) 

When  Tj  reaches  r',  this  slip  system  slides  and 
generates  e'l,.  Slip  then  occurs  both  in  the  ay?  slip 
system  in  P  and  Q,  and  in  the  system  in  R.  The 
incremental  resolved  shear  stresses  are  then 


ATa^.=  -  C{m,  a/5;  n,  a/5)Ae''^„ 

-C{m,aP;p,  |j/)a4,. 

(11) 

-C{m,^ri;p, 

(12) 

(The  repetitions  of  the  subscript  n  and  p  denote  sum¬ 
mation.)  Based  on  the  method  described  above,  the 
extrusion  growth  of  a  single  fatigue  band  with  loading 
cycles  was  calculated  [5,  6,  16]  for  the  case  with  no 
secondary  slip,  corresponding  to  an  ordered  super¬ 
alloy,  and  that  with  second  slip,  corresponding  to  the 
disordered  condition. 

We  have  also  considered  the  case  where  the  most 
favorable  crystal  has  11  fatigue  bands  [5].  During 
loading,  the  excessive  shear  stress  is  relieved,  owing  to 
slip  in  the  different  slip  bands.  Slip  in  one  band  relieves 
the  shear  stress  not  only  in  this  band  but  also  in  other 
bands.  The  center  band  has  the  maximum  initial  shear 
stress  and  gives  the  maximum  plastic  shear  strain  at  the 
free  surface.  This  shear  strain  is  considered  to  be  a 
measure  of  the  fatigue  damage.  The  excessive  shear 
stress  is  defined  as 

r^=T'  +  r^-T^  (13) 

The  variation  of  the  surface  plastic  strain  with 
loading  cycles  for  different  values  of  excessive  shear 
stress  for  both  the  ordered  and  disordered  inter- 
metallic  compounds  was  calculated  as  shown  in  ref.  16. 
From  these  results,  we  pick  the  number  of  cycles  to 
yield  the  given  value  of  the  surface  plastic  strain  for 
different  excessive  shear  stresses.  This  surface  plastic 
strain  is  a  measure  of  the  extent  of  extrusion  or  intru¬ 
sion  and,  hence,  the  extent  of  crack  initiation.  These 
vs.  N  curves  for  both  the  ordered  and  disordered 
superalloys  for  stress-controlled  loadings  are  shown  in 
Fig.  6.  It  can  be  seen  that  this  high  cycle  fatigue  crack 
initiation  life  in  stress-controlled  loading  of  the  ordered 
intermetallic  compounds  is  longer  than  that  of  the 
disordered  case  [15].  This  explains  the  test  results  of 
Boettner  et  al.  [2j. 

The  rate  of  incremental  slip  in  P  and  Q  per  cycle 
depends  on  the  resolved  shear  stress  (r'+rO-  As 
loading  progresses,  extrusion  increases,  causing 


( r'  +  rO  to  decrease  and  hence  reducing  the  growth 
rate.  This  decreases  the  slope  of  e'„|,  at  the  free  surface 
with  N.  For  the  multiple  fatigue  bands,  slip  in  one  band 
relieves  the  shear  stress  not  only  in  this  band  but  also  in 
other  bands.  This  causes  the  e'^pvs.  N  curve  to  bend 
further  down. 

4.  ].  Strain-controlled  cyclic  loading 

Under  strain-controlled  fatigue,  the  incremental 
macroscopic  plastic  strain  A£"  per  cycle  is  constant. 
AE"  was  taken  to  be  proportional  to  the  incremental 
surface  shear  strain  Ac"„^  in  the  case  of  a  single  fatigue 
band  in  the  most  favorably  oriented  crystal  and  to  be 
proportional  to  the  sum  of  the  incremental  surface 
strains  of  all  bands  in  the  case  of  multiple  fatigue  bands. 
The  method  of  calculation  is  the  same  as  in  the  case  of 
stress-controlled  loadings  [5,  6,  16].  The  applied 
stress— also  called  the  fatigue  resistance— was  varied  to 
yield  the  given  incremental  palstic  strain.  Two  values  of 
(second  fatigue  band)  or  the  sum  of  Ac"  ^  values 
(multiple  fatigue  bands),  i.e.  one  of  0.010  and  the  other 
of  0.026,  were  analysed.  The  variation  of  this  stress  vs. 
loading  cycles  is  shown  in  Fig.  7.  It  can  be  seen  that  this 
stress  Vi.  loading  cycles  curve  of  the  ordered  case  is 
higher  than  that  of  the  disordered  case.  This  agrees 
with  the  experimental  results  shown  in  Fig.  2. 

4.2.  Combined  cyclic  axial  and  torsional  loadings 
The  above  analysis  of  the  effect  of  long-range  order 
on  fatigue  crack  initiation  under  cyclic  axial  loadings  is 
here  extended  to  combined  cyclic  axial  and  torsional 
loadings.  Let  us  consider  a  circular  shaft,  with  x  along 
the  axial,  r  along  the  radial  and  6  along  the  circum¬ 
ferential  direction.  Single-crystal  tests  [18]  have  shown 
that  extrusions  occurred  on  the  slip  plane  along  the 
most  highly  stressed  slip  direction  and  did  not  occur 
when  the  direction  was  parallel  to  the  free  surface.  In 
cyclic  torsion,  the  maximum  shear  stress  occurs  along 
the  circumferential  direction,  parallel  to  the  free 
surface,  on  a  plane  normal  to  the  shaft  axis.  Hence,  the 
extrusion  or  intrusion  process  will  not  occur  in  the  slip 
system  with  maximum  shear  stress  but  may  occur  in 


Fig.  8.  Combined  axial  and  shear  loading. 


T.  H.  Lin,  Q.  Chen  /  Micromechanics  of  stress-  and  strain-controlled  fatigue 


109 


some  crystal  with  a  slip  plane  normal  to  the  shaft  axis 
and  a  slip  direction  making  an  angle  /?  with  the  circular 
boundary.  Under  combined  axial  and  torsional  loading, 
the  maximum  combined  stresses  occur  in  the  outer 
layer  of  the  shaft  and  gives  two  principal  stresses  a,  and 
an,  as  shown  in  Fig.  8.  The  third  principal  stress  am  is 
along  the  radial  direction  and  is  zero.  Referring  to  the 
X,  rand  d  coordinates,  we  have 


^  ^  _  ''XX 

- 


(14) 


The  three  extreme  values  of  shear  stresses  are  then 


surface  plastic  slip  strain  e",  with  and  without  the 
secondary  slip  were  calculated  for  different  values  of  /3. 
It  was  found  that  this  surface  plastic  strain  reaches  a 
maximum  at  about  /S  =  25°.  At  this  ^  value,  the  surface 
strain  e"^  with  secondary  slip  is  17%  more  than  that 
with  no  secondary  slip.  The  variation  of  this  plastic 
strain  with  loading  cycles  in  stress-controlled  cycling  is 
shown  in  Fig.  11.  This  gives  the  increase  in  surface 
plastic  strain  under  cyclic  loading  of  the  disordered 
above  that  of  the  ordered  intermetallic  compounds. 

Under  combined  axial  and  torsional  loadings,  the 
crack  initation  can  be  caused  by  T2  or  T3,  which  has  the 


where  r2  and  T3  have  a  shear  direction  making  an  angle 
of  45°  with  the  free  surface,  and  give  the  same  inclina¬ 
tion  angle  as  specimens  under  cyclic  tension  and  com¬ 
pression.  Hence,  the  analysis  for  cyclic  tension  and 
compression  is  applicable  to  evaluate  the  fatigue  crack 
initiation  caused  by  T2  and  T3.  However,  T|  has  the 
maximum  shear  stress  direction  parallel  to  the  free 
surface.  The  extrusion  and  intrusion  process  will  not 
occur  in  this  maximum  shear  stress  direction.  Never¬ 
theless,  this  process  can  occur  in  some  crystal  with  a 
slip  direction  making  an  angle  ^  with  the  free  surface 
(Fig.  9)  as  shown  by  Cooley  and  Lin  [  19]. 

Let  us  consider  such  a  crystal  which  has  a  slip  plane 
with  a  normal  along  the  x,  direction  (Fig.  10).  Thus,  we 
have 

T,^=T,  cos/5  (18) 

Two  closely  located  thin  slices  P  and  Q  which  have  a 
slip  direction  that  makes  an  angle  /3  with  the  free 
surface  are  assumed  to  exist  in  a  polycrystal,  as  shown 
in  Fig.  9.  As  in  the  cases  of  cyclic  tension  and  compres¬ 
sion,  positive  slip  in  P  and  negative  slip  in  Q  increase 
monotonically  with  loading  cycles  to  initiate  an  extru¬ 
sion  or  intrusion.  The  displacement  induced  by  the 
plastic  strain  is  along  the  ^  direction.  The  com¬ 
ponent  of  the  displacement  normal  to  the  free  surface 
Ufl  =  tel/3  is  taken  as  a  measure  of  the  amount  of  extru¬ 
sion  or  intrusion  and  the  crack  initiation. 

The  radial  component  u,  =  sin  P=te"p^  sin  p.  This 
Mr  at  the  free  surface  causes  an  extrusion  and  a  stretch 
in  R,  and  activates  the  secondary  slip.  This  secondary 
slip  in  R  will  then  increase  the  slip  in  P  and  Q.  The 


Fig.  9.  Cross-section  of  the  active  glide  plane.  The  slip  direction 
P  forms  an  angle  ^  with  the  free  surface  [19]. 


r 


Fig.  1 0.  Surface  crystal  with  two  slip  lines  P  and  Q  [  1 9). 


Honber  of  Cycles 

Fig.  11.  Surface  plastic  strain  vs.  number  of  cycles  of  an  inter¬ 
metallic  compound  under  a  cyclic  combined  axal  and  shear 
loadings. 
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same  extrusion  characteristic  as  that  under  cyclic 
tension  and  compression,  as  shown  in  the  first  part  of 
the  paper,  or  that  caused  by  r, ,  as  shown  in  this  section. 
The  fatigue  crack  initiation  is  controlled  by  the  highest 
extrusion  or  intrusion  caused  by  these  maximum  shear 
stresses  T|,  T2  and  T3. 


5.  Conclusion 

This  paper  gives  a  method  to  analyse  the  effect  of 
the  long-range  order  of  intermetallic  compounds  on 
high  cycle  fatigue  crack  initiation  under  axial  stress- 
controlled  and  strain-controlled  loadings.  The  effect  of 
these  intermetallic  alloys  under  combined  cyclic  axial 
and  torsional  loadings  subject  to  stress-controlled 
cycling  is  also  shown.  This  seems  to  explain  the  experi¬ 
mental  results. 
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HIGH-CYCLE  FATIGUE  CRACK  INITIATION  OF 
INTERMETALLIC  COMPOUNDS 
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EFFECT  OF  ELASTIC  ANISOTROPY  ON  FATIGUE  BAND 
IN  POLYCRYSTALLINE  SOLIDS 
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ABSTRACT 

Extrusions  and  intrusions  in  fatigue  specimens  are  favorable  sites  for  fatigue  cracks.  These  extnisions 
and  intrusions  have  been  considered  to  be  composed  of  three  thin  slices  P,  Q,  and  /{  in  the  surface 
crystal  of  a  polycrystal.  These  thin  slices  arc  taken  to  be  in  the  most  favorably-oriented  crystal  at  the 
free  surface.  Previous  siu(^  (Lin  and  Ito,  1969,  Lin,  1992)  concerns  mainly  the  fatigue  bands  in 
aluminum  and  its  alloys.  The  anisotropy  of  the  elastic  constants  of  these  alloys  is  insignificant  and 
accordingly  neglected.  Then  the  polyciystal  is  considered  as  an  elastically  isotropic  and  homogeneous 
medium.  However,  the  anisotropy  of  the  elastic  constants  of  some  other  metallic  crystals,  such  as 
titanium  and  some  intermetallic  alloys,  may  not  be  negligible.  In  the  present  paper,  the  effect  of 
anisotropy  is  analyzed. 

The  elastic  anisotropy  of  the  surface  crystal  is  considered.  Since  the  metal  exterior  to  this  surface 
crystal  is  composed  of  crystals  with  random  orientations,  it  is  assumed  that  the  surface  crystal  is 
embedded  in  an  elastically  isotropic  and  homogeneous  semi-infinite  medium.  Thus  it  is  reduced  to  an 
inclusion  problem.  This  heterogeneous  solid  is  transformed  into  a  homogeneous  one  by  using  Eshelby*s 
Equivalent  Inclusion  Method  (Eshelby,  1957,  1961,  Mura,  1982).  In  this  method,  an  extra  set  of 
eigens trains,  which  are  the  same  as  the  inelastic  strains,  is  introduced.  The  cigenstrains  explicitly 
depend  on  the  applied  stresses  and  plastic  strains  and  thus  can  be  implicitly  incorporated  in  the 
formulation  of  solutions.  The  remaining  calculation  procedures  arc  the  same  as  those  reviewed  by  Lin 
(1992).  The  anisotropic  elastic  constants  of  Ni^Al  monocrystalline  intermetallic  compoimd  (Yang, 
1985)  arc  used.  The  isotropic  elastic  constants  of  the  surrounding  metal  arc  taken  from  the  conventional 
average  method.  The  plastic  strain  distributions  and  the  cumulative  surface  plastic  strain  in  the  fatigue 
band  versus  the  ntimber  of  loading  cycles  are  calculated,  and  the  effect  of  elastic  anisotropy  on  the 
formation  of  extrusions  and  intrusions  is  examined  in  this  paper. 
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INTRODUCTION 

The  process  of  fatigue  failure  in  ductile  metals  can  be  classified  into  the  following  three  stages:  the 
track  initiation,  the  growth  and  coalescence  of  microcracks  to  form  some  macrocracks  and  the  propaga¬ 
tion  of  dominant  cracks  leading  to  the  final  catastrophic  failure.  In  the  hi^-cycle  fatigue,  most  of  fa¬ 
tigue  life  of  materials  is  spent  on  crack  initiation  which  starts  with  fatigue  bands.  Fatigue  bands  are  the 
favorable  sites  for  the  formation  of  fatigue  microcracks. 

Fatigue  bands  consist  of  intensely  localized  plastic  strains  and  have  been  observed  in  both  monocrys¬ 
talline  and  polycrystalline  metals.  Plastic  strains  are  mainly  induced  ly  the  displacement  trf 
dblocations.  Under  q^clic  loadings,  the  localized  plastic  d^ormation  accumulates  with  the  increase  rf 
the  number  of  loading  cycles  and  causes  extrusions  and  intrusions  on  the  surface  trf  materials.  The 
extent  of  extrusions  and  intrusions  is  here  taken  as  a  measure  of  fatigue  damage. 

A  ratchet  mechanism  has  been  proposed  to  model  the  fatigue  band  in  polyaystals  (Lin  and  Ito,  1969, 
Lin,  1992).  It  is  assumed  that  a  fatigue  band  composed  of  three  thin  slices  P,QmAR  is  embedded  near 
the  free  surface  in  an  elastically  isotropic  and  homogeneous  semi-infinite  medium.  Since  this  model  is 
concerned  with  the  fatigue  cradc  initiation  of  aluminum  and  its  alloys,  the  elastic  anisotropy  of  their 
constituent  monocrystals  is  insignificant  and  accordingly  neglected.  For  some  other  polycrystalline 
metals,  such  as  titanium  and  some  intermetallic  alloys,  however,  the  elastic  anisotropy  of  crystals  may 
not  be  negligible.  Since  the  crystal  orientations  can  vary,  the  local  anisotropy  and  inhomogeneity  of 
polycrystals  will  affect  the  formation  of  the  fatigue  band. 

In  this  paper,  the  elastic  anisotropy  of  the  surface  crystal  is  considered.  Since  the  metal  exterior  to  the 
surface  crystal  is  composed  of  crystals  with  random  orientations,  the  surface  crystal  is  taken  as  being 
embedded  in  an  elastically  isotropic  and  homogeneous  semi-infinite  medium.  Thus  this  analysis  is 
reduced  to  an  inclusion  problem.  This  heterogeneous  solid  is  transformed  mto  a  homogeneous  one  by 
using  the  Eshelby’s  Equivalent  Inclusion  Method  (Eshelby,  1957, 1961,  Mura,  1982).  hr  this  method,  an 
extra  set  of  eigentrains,  whidi  are  the  same  as  the  inelastic  strains,  is  introduced.  These  eigenstrains 
can  be  explicitly  expressed  in  terms  of  applied  stresses  and  plastic  strains  as  well.  After  this 
transformation,  the  procedures  of  analysis  are  the  same  as  those  reviewed  by  Lin  (1992).  fri  the 

following  numerical  calculations,  the  anisotropic  elastic  constants  trf  Ni^Al  monocrystalline  interme¬ 
tallic  compound  (Yang,  1985)  are  used.  The  isotropic  elastic  constants  of  the  surrounding  metal  are 
taken  from  the  conventional  average  method.  The  plastic  strain  distributions  and  the  cumulative  surface 
plastic  strains  in  the  fatigue  band  versus  the  number  of  loading  cycles  are  calculated,  and  the  effect  rf 
elastic  anisotropy  on  the  formation  of  extrusions  and  intrusions  is  examined. 

CALCULATION  OF  RESIDUAL  STRESS  HELD 

The  residual  stress  field  due  to  a  distribution  of  inelastic  strains  can  be  calculated  by  the  analogy 
between  inelastic  strains  and  applied  forces  (Lm,  1969).  It  has  been  shown  that  the  equivalent  body 

force  per  unit  volume  along  the  x,-  axis  due  to  inelastic  strains  e", ^  is 

Fi  =  -C^^€'u.j  (1) 

where  C,^  is  the  elastic  modulus.  The  summation  convention  is  used  here  and  the  subscript  after  a 

comma  denotes  the  differentiation  with  respect  to  the  coordinate  variable.  For  isotropic  elastic  solids, 
the  elastic  constants  are 

+ p.(5jj5yj  +  5,-^  ^) 
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where  X  and  are  Lame’s  constants  and  5-  is  the  Kronecker  delta.  The  equivalent  surface  force  per  unit 
area  along  the  Xf-  axis  has  been  shown  as 

5/  =  Cj,A.,e''*,Vy  (2) 

where  v^.  is  the  direction  cosine  of  the  normal  to  the  area.  The  fictitious  stresses  a-fx)  at  point  x 
caused  by  equivalent  forces  in  domain  Q  can  be  written  by  using  the  Green's  function  as 

oJ(x)  =  J  t|(x.x’  )F*(x’  )rfQ+ jigCxx-  )Sjt(x'  )rfr 

a  r 

t  t  r  . 

^WmnJ  Q./»nnJ 

a  r 

where  t|(x,x’  )  represents  the  stress  component  0,^-  at  point  x  caused  by  a  unit  force  acting  at  point  x’ 
along  the  -axis.  P  is  the  boundary  of  domain  £2.  The  residual  stress  field  is 

aJ(x)=of(x)-Ci;*,8"*,(x)  (4) 

To  calculate  the  residual  stress  field  numerically,  the  discretization  of  a  continuous  problem  is 
needed.  In  the  following,  the  donnain  £2  in  an  extended  solid  is  divided  into  a  number  of  subdomains 

In  each  of  them,  it  is  assumed  that  the  inelastic  strains  e*,^-  are  uniform.  After  applying  eqns  (3)  and 

(4)to  n„,  the  residual  stress  field  due  to  £"„,„(x')(x’€not)  alone  is 

j nin(^a) 

r« 

where  the  function  //(x,Q„)  is  defined  as 


H(x,n  )= 


Taking  the  average  of  crffx)  in  subdomain  flp  as 


•  1,  xeQ„ 

-  0,  xsQo^ 


JoJ(x)f/Q 

and  letting  o-j(Qp)  represent  the  residual  stresses  in  Hp  due  to  in  £2n,  then  we  obtain 

~  ^ijntn  (f^p>^a in/i  (^o ) 


(5) 
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'^'lere  are  called  the  residual  stress  influence  coefficients,  which  represent  the  residual 

stresses  in  subdomain  caused  by  unit  inelastic  strain  components  in  subdomain  n  . 


GENERALIZATION  OF  EQUIVALENT  INCLUSION  METHOD 

valent ff"'  by  Eshelby  (1957,  1961).  which  is  now  called  the  Eshelby’s  Equi- 

t  Inclusion  Method,  can  be  used  to  transform  an  inhomogeneous  problem  into  a  homogeneous  one 
by  intrc^ucing  a  distribution  of  eigenstrains  (Mura,  1982).  Originally,  this  method  was  appL  to  solve 

problems  with  an  ellipsoidal  inhomogeneity.  In  this  study,  this  method  is  generalized  to  solve 
inhomogeneous  problems  of  various  shapes. 

In  the  following,  the  matrix  and  vector  notation  will  be  used  instead  of  tensor  notation  for  con- 

r  T?'*  boldface 

uppercase  Roman  letters.  The  Eshelby's  equivalent  equation  is 

C*(x)(e“  +  E(x)-e/'(x))  =  C(e‘‘  +  £(x)-e*(x)-e/’(x))  .  xeQ  (6) 

where  Q  is  an  inclusion  or  a  subdomain  with  the  elastic  constant  C*(x)  embedded  in  an  extended  body 

Z)  with  the  isotropic  elastic  constant  C.  e(x)  is  the  disturbance  strain,  8*(x)  the  eigenstrain,  eZ’(x) 

the  plastic  strain  and  e“  uniform  strain  from  which  <r*  =  Ce“  where  a"  is  the  applied  far-field  stress. 

Eqn(6)  shows  that  the  eigenstrain  £*(x)  plays  the  same  role  as  the  plastic  strain  £/’(x)  does  in  the 
homogeneous  materials. 

When  a  continuous  problem  is  discretized  for  the  purpose  of  numerical  analysis,  the  domain  of  Q  is 

further  divided  into  N  subdomains  Q.,(\<  i<N)  called  elements  in  which  both  of  the  eigenstrain  e*and 

plastic  strain  e'’  are  assumed  to  be  uniform.  The  stresses  in  element  Q,  caused  by  the  unit  eigenstrain 

e  in  element  Qj  are  related  by  the  residual  stress  influence  coefficients,  the  same  as  caused  by  the 

plastic  strain  e^.  Thus  the  stress  of  in  Q,  induced  by  both  of  the  eigenstrains  and  plastic  strains  in  all 
the  elements  concerned  is 


where 


®;  =  IG,e;+XG„£f  , 

]  k 


1  SiSN 


(7) 


^22* 

^33'  ^23’ 

^!13» 

J 

^22- 

pP  pP 

^33’  ^23’ 

erj; 

in  Cl.  and  are  6x6  matrices  of  the  residual  stress  influence  coefficients,  which  are  contracted  from 

Gkimn(«i.Gj).  Notation  {-f  represents  the  transpose.  Since  the  generalized  plane  strain  problem  is 
dealt  with  in  this  study,  the  following  equation  also  holds 

of  =  Ce,.-C(e;  +  8f) 


where 


(8) 
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£.^=0 


ei={eil.  £22*  £23»  £13*  £12},^  with 
From  eqns(7)  and  (8),  solving  for  disturbance  strains  e,-  in  £2;  in  terms  of  ej  and  we  obtain 


where 


£;  =  C"' 


C(e;+ef)+IG^£}  +  IG,.tE,^ 


lSi<N 


(9) 


s 

d 

n 

e 

I 


C-'=1 


1  -V  -V  0  0  O' 

-V  1  -V  0  0  0 

000  1+v  0  0 

0000  1+v  0 

0  0  0  0  0  1  +  v 


in  which  £  is  the  Young’s  modulus  and  v  the  Poisson's  ratio. 

The  left-hand  side  of  eqn(6)  can  also  be  written  for  the  generalized  plane  strain  problems  as 


<y,=0:+C*8,-C*ef 


(10) 


where  C-"  is  a  6x6  anisotropic  elastic  stiffness  matrix  and  C*a  6x5  matrix  formed  from  matrix  C*  by 
taking  the  third  column  away.  And  a*  =  C*£'’  in  which 

®  f  O  o  o  1 

£  ==  \£ii,  £22*  £12/  with  £33*0 

Substituting  eqn(9)  into  eqn(lO)  and  eqn(7),  eqn(6)  becomes 


where 


Te*  +  =  q;  +  (C*  -T)ef  - 

j  k 


1  <i  <N 


(11) 
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T  =  C*C"'C 

o  * 

qi=0  -O, 

(y°  =  C£'’ 

and  I  is  a  6x6  unitmatrix. 

It  can  be  seen  that  eigenstrains  £*  depend  on  the  applied  load  a°  and  the  distribution  of  plastic 
strains  E]^  as  well. Solving  eqn(l  1)  for  t*  in  allelements  Q|,  we  have 

e*  =P/+XR;jtelt  >  1^‘^N  (12) 

k 

where  p,-  is  a  vector  in  terms  of  o'  in  Q, .  The  total  stress  in  Q; ,  obtained  by  substituting  eqn(l  2)  into 
eqn(7),  is 

0,=a“+o[ 

=o”+IGjPj+I(G«+XGjR;,)ef 
J  k  j 

In  the  above  equation,  the  term  of  £GjjPj  is  a  modification  to  the  effect  of  the  applied  load  o°  and 

X^ijRji;  the  one  to  the  matrix  Gj^  of  the  residual  stress  influence  coefficients  due  to  the  inhomo¬ 
geneity. 

If  there  is  only  the  stress  component  o^  applied,  then  it  is  denoted  in  eqn(13)  that 

P/cy22 

j 

and 

^ik  =  +  X  Gy-  R  Jic 

j 

Eqn(l  3)  becomes 

k 

For  slip  system  ap  in  element  Qj ,  the  following  relations  hold 

ef  =  “apePpi 

'^api "  ‘•op 
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where  e^pj  and  are  the  plastic  shear  strain  and  resolved  shear  stress,  respectively,  and 

=  «*22.  «»33>  "^23’  «Jl3> 

"ap~{"n’  '*22>  '*33>  ”23'  ”l3,  ”12]qP 

in  which  mjj  and  njj  are  composed  of  direction  cosines  of  slip  system  ap.  The  total  resolved  shear  stress 
in  the  slip  system  is 

^op/  =  tap,+4/+^i, 

where 

=  "oip  Pi^h  ~ 

k  ^Ti 

in  which 

C(0(P,Q;  ;^T|,Qj(. 

and  the  sumnialion  with  regard  to  slip  system  includes  all  the  slip  systems  in  i^-th  element  in 
question. 

•^a|if  are  the  modified  Schmid  factors  and  the  residual  stress  influence  coeffi¬ 

cients  ,  respectively,  after  transforming  the  inhomogeneous  problem  into  a  homogeneous  one. 

PHYSICAL  MODEL  AND  NUMERICAL  CALCULATION 

A  polycrystal  model  has  been  established  to  study  the  effect  of  crystal  anisotropy  and  inhomogeneity 
on  the  fatigue  crack  initiation  as  shown  in  Fig.  1.  An  elastically  anisotropic  crystal  is  embedded  near  the 
free  surface  in  a  semi-infinite  elastic  medium.  Since  the  orientations  of  surrounding  crystals  are  random, 
it  is  assumed  that  the  metal  exterior  to  the  surface  crystal  is  isotropic  and  homogeneous.  A  fatigue  band 
composed  of  three  thin  slices  Q  and  R  is  in  this  surface  crystal,  which  is  taken  as  to  be  the  most 
favorably-oriented  with  the  primary  slip  system  being  45°  inclining  to  the  load  axis. 

The  elastic  constants  of  a  nickel-base  superalloy  Ni^Al  single  crystal  with  cubic  symmetry  (Yang, 

1 985)  are  used  in  the  numerical  calculation.  The  elastic  compliance  moduli  in  10'“/CPa  with  respect  to 
the  crystal  axes  [001],  [010]  and  [100]  in  the  same  notation  as  Yang's  are 

5l,  =  1.01 

5,2  =-0.393 

544=0.848 

The  shear  modulus  for  the  isotropic  and  homogeneous  material  is  75  GPa,  which  is  obtained  by  taking 
the  average  stress  over  a  number  of  orientations  of  the  single  crystal  corresponding  to  a  given  uniform 
strain,  and  the  Poisson's  ratio  is  0.3. 

In  order  to  fit  the  geometry  of  the  surface  crystal  and  the  fatigue  band,  a  constant  strain  trapezoidal 
element  in  general  and  a  linear  strain  element  for  the  fatigue  band  in  particular  arc  developed  to 
calculate  the  residual  stress  influence  coefficients  for  the  generalized  plane  strain  problem.  The  area  of 
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the  surface  crystal  is  divided  into  50  constant  strain  elements  and  each  slice  of  P,  Q  and  R  into  12  linear 
strain  elements. 

With  the  method  described  previously,  the  modified  Schmid  factor  varying  along  slice  P  is  shown  in 
Fig.  2.  Compared  with  the  case  without  considering  the  crystal  anisotropy  where  the  Schmid  factor  is 
uniformly  0.5,  the  anisotropy  and  inhomogeneity  lower  this  value.  Since  the  slip  theory  is  used  in  this 
study,  the  slip  system  in  which  the  total  resolved  shear  stress  reaches  the  critical  shear  stress  may 
slide  first.  Therefore,  the  sliding  starts  inside  the  band,  where  the  modified  Schmid  factor  takes  the 
maximum  value.  The  development  of  fatigue  band  depends  upon  the  irreversibility  of  the  plastic  shear 
deformations  in  P  and  Q.  This  irreversibility  of  sliding  in  the  fatigue  band  is  caused  by  an  initial  stress 
field  T  ,  which  is  induced  by  an  initial  strain  field  (Lin  and  Lin,  1983).  In  the  numerical  calculations, 
the  initial  stress  x'  is  0.5  MPa.  The  critical  shear  stress  x^^  is  taken  as  lOOMPa  and  the  applied  loading 
magnitude  is  AZOMPa.  Since  the  sliding  outside  the  fatigue  band  is  reversible  during  the  cyclic 
loading  process,  the  plastic  deformation  concentrates  in  the  fatigue  band.  Fig.  3  shows  the  extrusion 
growth  of  the  fatigue  band  with  cycles  of  loading.  The  plastic  shear  strain  distributions  in  P  and  Q  are 
plotted  at  various  loading  cycles  in  Fig.  4. 


CONCLUDING  REMARKS 

A  generalized  Eshelby's  Equivalent  Inclusion  Method  is  developed  to  study  the  effect  of  the  elastic 
anisotropy  of  crystals  on  the  crack  initiation  of  polycrystalline  metals.  Since  this  is  a  local  effect,  the 
results  depend  on  the  elastic  properties  of  both  the  surface  crystal  and  the  sunounding  metal.  In  this 
study,  the  elastic  constants  for  the  surrounding  metal  are  of  the  upper  bound  by  nature.  The  selection  of 
these  elastic  constants  should  reflect  the  average  effect  of  the  surrounding  material  on  the  surface 
crystal.  A  more  realistic  model  which  incorporates  the  surrounding  crystals  with  different  orientations  is 
being  developed  to  study  this  effect. 
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Fig.  3  Cumulative  Plastic  Shear  Strain  at  the  Free  Surface  in  P  vs.  the  Number  of  Cycles 


Fig.  4  Cumulative  Plastic  Strain  Distributions  in  P  and  Q  at  Cycles  of  Loading 
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Fatigue  bands  have  been  observed  in  both  monocrystalline  and  polycrystalline  metals. 
E^rusions  and  intrusions  at  the  free  surface  of  fatigued  specimens  are  favorable 
sites  for  fatigue  crack  nucleation.  Previous  studies  {Lin  and  ho,  1969;  Lin,  1992) 
mainly  concerned  the  fatigue  crack  initiation  in  aluminum  and  its  alloys.  The  elastic 
anisotropy  of  individual  crystals  of  these  metals  is  insignificant  and  was  accordingly 
neglected.  However,  the  anisotropy  of  the  elastic  constants  of  some  other  metallic 
crystals,  such  as  titanium  and  some  intermetallic  compounds,  is  not  negligible.  In 
this  paper,  the  effect  of  crystal  anisotropy  is  considered  by  using  Eshelby*s  equivalent 
inclusion  method.  Hie  polycrystal  analyzed  is  NisAl  intermetallic  compound.  The 
plastic  shear  strain  distributions  and  the  cumulative  surface  plastic  strain  in  the 
fatigue  band  versus  the  number  of  loading  cycles  were  calculated,  and  the  effect  of 
crystal  anisotropy  on  the  growth  of  the  extrusions  was  examined. 


1  Introduction 

The  process  of  fatigue  failure  in  ductile  metals  is  divided 
into  the  following  three  stages:  ( 1 )  the  crack  initiation,  (2)  the 
growth  and  coalescence  of  ihicrocracks  to  form  some  mac¬ 
rocracks,  and  (3)  the  propagation  of  dominant  cracks  leading 
to  the  final  catastrophic  failure.  In  high-cycle  fatigue,  fatigue 
crack  initiation  predominates  fatigue  life.  Fatigue  bands  are  the 
favorable  sites  for  the  nucleation  of  fatigue  microcracks. 

Fatigue  bands  consist  of  highly  localized  plastic  strains  and 
have  been  observed  in  both  single  crystals  and  polycrystals. 
Plastic  strains  ^e  mainly  due  to  the  movement  of  dislocations. 
Under  the  cyclic  loading,  the  localized  plastic  deformation  in¬ 
creases  with  the  number  of  loading  cycles  and  causes  extrusions 
and  intrusions  at  the  surface  of  metals.  The  height  of  extru- 
sions  or  the  depth  of  intrusions  is  here  taken  as  a  measure  of 
fatigue  damage. 

A  theory  of  micromechanics  of  fatigue  crack  initiation  has 
been  proposed  to  model  the  gating  mechanism  of  fatigue  bands 
in  polycrystals  (Lin  Ito,  1969;  Lin,  1992).  A  fatigue  band 
composed  of  three  thin  slices  P,  Q,  and  P  in  a  most  favorably 
oriented  surface  crystal  was  considered.  Under  cyclic  tension 
and  compression  loading,  the  alternate  monotonic  build-up  of 
plastic  shear  straiiis  in  both  P  and  Q  causes  the  growth  of 
extrusions  or  intrusions  at  the  free  surface  of  polycrystals.  Since 
this  model  is  concerned  mainly  with  the  fatigue  crack  initiation 
of  aluminum  and  its  alloys,  the  elastic  anisotropy  of  their  indi¬ 
vidual  crystals  is  insignificant  and  accordingly  neglected.  Hence 
the  surface  crystal  is  assumed  to  be  embedded  in  an  elastically 
isotropic  and  homogeneous  semi-infinite  medium.  For  some 
other  polycrystalline  metals,  such  as  titanium  and  some  inter- 
metalic  superalloys,  however,  the  elastic  anisotropy  of  crystals 
is  not  negligible.  For  such  polycrystals,  the  local  microstress 
field  in  a  crystal  depends  not  only  on  the  applied  stresses  but 
also  on  the  orientations  of  its  surrounding  crystals.  Hence  the 
elastic  anisotropy  of  crystals  will  affect  the  fatigue  crack  initia¬ 
tion  in  these  metals. 

In  this  paper,  the  elastic  anisotropy  of  individual  grains  in  a 
polycrystal  is  considered.  A  fatigue  band  is  assumed  to  occur 
in  the  most  favorably  oriented  surface  crystal.  This  crystal  is 
surroimded  by  the  crystals  of  various  given  orientations.  The 
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crystals  outside  the  surrounding  crystals  are  randomly  oriented 
and  assumed  to  be  isotropic  and  homogeneous.  Hence  this  sur¬ 
face  crystal  and  the  surrounding  crystals  of  given  orientations 
are  considered  to  be  embedded  in  a  semi-infinite  isotropic  and 
homogeneous  elastic  medium.  Thus  it  becomes  an  inclusion 
problem.  This  heterogeneous  solid  is  transformed  into  a  homo¬ 
geneous  one  by  using  Eshelby’s  equivalent  inclusion  method 
(Eshelby,  1957,  1961;  Mura,  1982).  In  this  method,  an  extra 
set  of  eigenstrains  is  introduced.  These  eigenstrains  can  be  ex¬ 
plicitly  expressed  in  terms  of  the  applied  stresses  and  plastic 
strains.  After  this  transformation,  the  procedures  for  this  micro¬ 
mechanical  analysis  are  the  same  as  those  given  by  Lin  (1992) 
for  elastically  isotropic  polycrystals.  In  the  following  numerical 
calculation,  the  anisotropic  elastic  constants  of  NisAl  monocrys- 
^line  intermetallic  compound  (Yang,  1985)  were  used.  The 
isotropic  elastic  constants  of  the  surrounding  metal  were  taken 
ftom  the  experimental  data  (Stoloff,  1989).  The  plastic  strain 
dis^butions  and  the  cumulative  surface  plastic  strains  in  the 
fatigue  band  versus  the  number  of  loading  cycles  were  calcu¬ 
lated,  and  the  effect  of  elastic  anisotropy  on  the  formation  of 
extrusions  and  intrusions  was  examined. 

2  Gating  Mechanism  Provided  by  Microstress  Field 

Sin^e  crystal  tests  show  that  slip  occurs  on  certain  crystallo¬ 
graphic  planes  along  certain  directions.  This  slip  depends  on 
the  resolved  shear  stress  in  the  slip  system  and  is  independent 
of  the  normal  stress  on  the  sliding  plane.  The  dependence  of 
slip  on  the  resolved  shear  stress  is  known  as  the  Schmid  law. 
It  has  been  shown  that  the  Schmid  law  also  holds  when  crystal¬ 
line  solids  are  subjected  to  cyclic  loading  (Parker,  1961).  For 
a  face-centered  cubic  crystal  (FCC),  there  are  four  slip  planes, 
on  each  of  which,  there  are  three  slip  systems.  These  twelve 
slip  systems  are  shown  in  Fig.  1. 

All  met^s  have  initial  defects,  hence,  have  initial  stresses, 
^ese  initial  stresses  play  an  important  role  in  fatigue  crack 
initiation.  Lin  and  Ito  ( 1969 )  first  point  out  that  the  initial  stress 
field  favorable  to  extrusions  and  intrusions  can  be  caused  by  a 
row  of  dislocation  dipoles.  Later,  these  dislocation  dipoles  have 
been  observed  in  the  ladder  structure  in  a  persistent  slip  band 
(PSB)  (i^tonopoulos,  1976;  Mughrabi,  1982).  These  dipoles 
induce  initial  resolved  shear  stresses  at  the  interface  between 
the  PSB  and  the  matrix. 

Based  on  experimental  observations,  a  xnicromechanic  theory 
of  fatigue  crack  initiation  was  proposed  by  Lin  and  Ito  ( 1969). 
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Slip  Systems  and  Their  Directions 
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Fig.  1  Crystallographic  slip  planes  and  directions  of  FCC  crystals 


In  this  model,  a  most-favorably  oriented  surface  crystal  con¬ 
taining  a  fatigue  band  composed  of  three  thin  slices  P,  Q,  and 
R  in  a  semi-infinite  polycrystal  shown  in  Fig.  2  is  considered. 
The  slip  plane  and  slip  direction  of  the  primary  slip  system  a/3 
form  an  angle  of  45  deg  with  the  loading  axis.  Lin  and  Lin 
(1983)  have  calculated  a  uniform  initial  resolved  shear  stress 
field  positive  in  P  and  negative  in  Q»  caused  by  a  distribution 
of  plastic  strains.  This  initial  stress  field  satisfies  the  conditions 
of  both  compatibility  and  equilibrium.  During  the  cyclic  load¬ 
ing,  plastic  strains  occur  in  P,  Q  smd  R  and  cause  a  residual 
stress  field.  The  total  resolved  shear  stress  t  in  a  slip  system  is 
the  sum  of  three  stresses  as 


where  is  the  resolved  shear  stress  caused  by  the  applied  load 
and  the  residual  resolved  shear  stress  due  to  plastic  strains. 

A  tensile  loading  cr  22  on  the  polycrystal  produces  a  positive 
in  the  entire  surface  crystal.  Since  rMn  F  is  positive,  the 
total  resolved  shear  stress  r  in  F,  which  is  +  r',  reaches 
the  critical  shear  stress  first;  and  hence,  according  to  Schmid 
law,  P  slides.  The  plastic  shear  strain  caused  by  this  slip  in  P 
produces  a  residual  stress  field  which  has  negative  resolved 
shear  stress  in  Q  and  makes  Q  to  slide  more  readily  in  the 
reverse  loading.  The  negative  slip  in  Q  induces  in  turn  a  positive 
residual  resolved  shear  stress  in  P  making  P  more  readily 
to  slide  under  the  next  tensile  loading.  This  process  is  repeated 
during  the  cyclic  loading  for  every  cycle  thus  providing  a  natural 
gating  mechanism  for  a  monotonic  buildup  of  local  resolved 
plastic  shear  strain  eS/3  in  both  P  and  Q,  pushing  R  out  of  the 
free  surface  and  promoting  an  extrusion.  The  elongation  of  R 
may  activate  a  secondary  slip  system  to  slide.  The  effect  of  slip 
in  the  secondary  slip  system  enhances  the  extent  of  extrusion 
(Lin,  et  al.,  1989).  Interchanging  the  signs  of  the  initial  stresses 
in  P  and  Q  initiates  an  intrusion  instead  of  an  extrusion.  This 
micromechamcal  theory  on  fatigue  crack  initiation  has  over¬ 
whelming  metallurgical  experimental  evidences  (Lin,  1977). 

3  Calculation  of  Residual  Stress  Field 

3.1  Residual  Stress  Field.  The  residual  stress  field  due 
to  a  distribution  of  inelastic  strains  can  be  calculated  by  the 
analogy  between  inelastic  strains  and  applied  forces  (Lin, 
1969).  It  has  been  shown  that  the  equivalent  body  force  per 
unit  volume  along  the  Xi-axis  due  to  inelastic  strains  €,yis 

~  (1) 

where  Ciju  is  the  elastic  modulus.  The  summation  convention 
is  used  here  and  the  subscript  after  a  comma  denotes  the  differ¬ 
entiation  with  respect  to  the  coordinate  variable.  For  isotropic 
elastic  solids,  the  elastic  constants  are 


Qjia  -  \6ij6u  +  +  6ii6jk) 

where  \  and  fx  are  Lame’s  constants  and  5,^  is  the  Kronecker 
delta.  The  equivalent  surface  force  per  unit  area  along  the  xr 
axis  has  been  shown  as 


n  n  f  1 1  n 


I  n  n  n  n 

Fig.  2  Fatigue  band  in  the  most  favorably  oriented  surface  crystal 
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Si  “  Cijki€ul^j  (2) 

where  Uj  is  the  direction  cosine  of  the  normal  to  the  area.  The 
fictitious  stresses  o’J(x)  at  point  x  caused  by  equivalent  forces 
in  domain  Cl  can  be  written  by  using  the  Green’s  function  as 

x')Ftix')dn  +  r|(x.  x’)Ft(x')dr 

=  -Cjj™  I  tJ(x,  x')€L./(x')<ff2 

+  J^rJ(x.x')€L(x')Mr  (3) 

where  r  ”(x,  x')  represents  the  stress  component  atj  at  point  x 
induced  by  a  unit  force  acting  at  point  x'  along  the  xt-axis.  T 
is  the  boundary  of  domain  If  inelastic  strain  components  of 
^ill«(x)  are  constant  in  fl,  eqn  (3)  is  reduced  to 

Tt(x,  x')i/,<ir|e;L.(n)  (4) 

where  £^(0)  is  the  constant  inelastic  strain  in  0.  The  residual 
stress  field  is 

(TijW  =  crf)(x)  -  Cijutui^)  (5) 
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Fig.  3  Euler  angles  and  representation  of  eiystal  orientations 


3.2  Rei^dual  Stress  bifluetice  CoefScieiiL  To  calculate 
the  residual  stress  field  numerically,  the  discretization  of  a  con¬ 
tinuous  problem  is  needed.  In  the  following,  the  domain  fl  in 
an  extended  solid  is  divided  into  a  number  of  subdomains  or 
elements  fl,.  In  each  of  them,  it  is  assumed  that  the  inelastic 
strain  e,*  is  uniform.  After  applying  Eqs.  (4)  and  (5)  to  fi,,  the 
residual  stress  field  due  to  e  fl,)  is 


~  ^Jj.  —  /f(x,  (6) 

where  the  function  H{x,  Cl,)  is  defined  as 


^f(x.  Cl,) 


_  fl,  X 
"to,  X 


e  n, 

« fl. 


Generally,  residual  stress  field  ajfx)  varies  from  one  point  to 
another.  The  volume  average  of  t7j(x)  ovct  subdomain  Clp  is 
taken  to  represent  the  residual  stresses  in  element  fl,.  Letting 
o'jiClf)  represent  the  residual  stresses  in  fl,  due  to  in  Cl,, 
we  have 


0»22 


1  )  »  J  1 


I  1  1  )  J 

0«22 

Fig.  4  Polycrystal  model  of  fatigue  crack  initiation 


(7) 

where  Gij^(Qpy  fi,)  are  called  the  residual  stress  influence  coef- 
flcients,  which  represent  the  residual  stresses  in  subdomain  Qp 
caused  by  unit  inelastic  strain  components  in  subdomain 
For  the  resolved  residual  shear  stress  in  slip  system 
in  due  to  the  resolved  plastic  shear  strain  ef,,  in  slip  system 
^77  in  n,,  Eq.  (7)  can  be  rewritten  by  the  tensor  transformation 
as 


rl0{np)  =  G(a/3,  ^t?,  (g) 

where  G(aP,  Qp ;  represents  the  residual  resolved  shear 

stress  in  slip  system  a0  in  flp  caused  by  a  unit  resolved  plastic 
shear  strain  in  slip  system  in  0^.  The  total  resolved  shear 
stress  Ta/j  in  slip  system  a0  in  f2p  is 

T.^Cl,)  =  T'.e  +  Tia  +  'Z1  Giotp,  Cl,-,  ^rj,  Cl,)e%, 

(9) 

The  summation  with  regard  to  slip  system  $77  includes  all  the 
slip  systems  in  fl,. 

33  Generalized  Plane  Strain  Problem*  The  thickness 
of  fatigue  bands  is  much  smaller  than  the  length  (dimension 
along  the  a:3-axis)  as  observed  at  the  free  surface  of  metals. 
Therefore,  the  problem  is  considered  under  the  plane  strain 
deformation.  The  plastic  strain  due  to  slip  in  the  secondary  slip 
system  is  not  confined  to  the  XiXz  plane.  This  causes  an  equiva¬ 
lent  force  component  acting  along  the  Xs-direction.  The  exis¬ 
tence  of  this  ^3  requires  a  modification  of  ihe  plane  strain  prob¬ 
lem.  Hence,  the  generalized  plane  strain  problem  is  considered 
in  the  present  study.  A  similar  problem  is  shown  in  the  analysis 
of  prismatic  anisotropic  bars  by  Lekhnitski  ( 1963) .  This  gener¬ 
alized  plane  strain  deformation  is  defined  as 

tt.  =  Ui(xuX2),  i  =  1,  2.  3 

where  w,-  is  the  displacement  component  along  Xf-axis.  This 
gives 

*  i  '^•'^1 

where  9  =  Mu  +  M2.2  ^  is  the  Poisson’s  ratio.  The  equilib- 
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Fig.  6  Modified  Schmid  factor  varying  along  slice  P 


Fig.  8  Plastic  shear  strain  distributions  in  P  at  100  cycles  and  2000 
cycles 


rium  equations  in  terms  of  displacement  components  u,-  can  be 
expressed  as 


1  5(9  F.  ^  _ 

1 — f  =  1,  2 
1-2/1  axi  fi 


(10) 


and 


and 


Tt3(x.x')  =  T|3(x.x')  =  0,  A:  =1,2  (12) 


ri3(x,  x')  = 


dx,  ' 


t|3(x,  x')  =^, 

dXi 


where 


(11) 


dx]  dxl 


and  Fi  is  the  body  force  component  Equations  (10)  and  (11) 
are  not  coupled  and  can  be  solved  separately. 

The  stress  comf^nents  r”(x,  x")  in  Eq.  (3)  for  the  semi- 
infinite  medium  with  the  generalized  plane  strain  deformation 
have  been  expressed  in  terms  of  the  Airy’s  stress  functions  (Lin 
and  Lin,  1974)  as 


r{,(x,  x')=^,  r|2(x,x')  =  ^, 
dxl  dxj 

Tt2(X,  x')  =  -  ^ 


dx,dx2 


,tU(x,  x')  = 


Tyix,  x')  =  0,  for  other  /  and  j  (13) 

where  the  Airy’s  stress  functions  are 
<^i(x,  X')  =  -(p  +  q){x2  -  x'lXe^  +  92) 

+  \qixt  -  •*;)  In  (Xi/Xj)  +  2pxix\(xt  +  x’iVXi 

«fe(X,  X')  =  (p  +  ^)(X2  -  X2)(^l  +  ^2) 

+  ^q(X2  -  Xl)  to  (X./Zj)  -  2pXiX\{X2  -  X’2)/X2 
<^,(x,  X')  =  -  ^  (to  X,  +  to  Xj) 


with 


^  47r(l 

- q=p{\-2v) 

- 

dx  =  arctan  ( 

\ 

\Xi  —  / 

02  =  arctan  ( 

\ 

Xi  +x\J  2  2 

Xi  =  (JC,  -  xD*  +  (X2  -  x^)^ 

X2  =  (Xi  +  x{)^  +  (.X2  -  X2)* 

By  using  Eqs.  (6),  (12)  and  (13),  the  residual  stresses  as 
well  as  the  residual  stress  influence  coefficients  due  to  the  in¬ 
elastic  strains  in  the  element  of  a  specific  shape  can  be  readily 
calculated. 

4  Generalization  of  Equivalent  Inclusion  Method 

Eshelby’s  equivalent  inclusion  method  (Eshelby,  1957, 
1961 )  is  here  used  to  transform  an  inhomogeneous  problem  into 
a  homogeneous  one  by  introducing  a  distribution  of  eigenstrains 
(Mura,  1982).  Originally,  this  method  was  applied  to  solve  the 
problems  with  an  ellipsoidal  inhomogeneity.  In  the  present 
study,  this  method  is  generalized  to  solve  inhomogeneous  prob¬ 
lems  of  arbitrary  shapes. 
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In  the  following,  the  matrix  and  vector  notation  will  be  used 
instead  of  tensor  notation  for  convenience.  It  is  easy  to  contract 
a  second-order  tensor  with  components  Gi^ta  into  a  matrix  with 
entries  Gij  by  setting,  for  instance,  Gu  =  Gmu  Gn  =  G,,22, 
Gi6  =  Ghi2  +  G1121  and  Ggi  =  G^n,  etc.  Vectors  are  denoted 
by  boldface  lowercase  Roman  or  Greek  letters  and  matrices  by 
boldface  uppercase  Roman  letters.  The  Eshelby’s  equivalent 
equation  is 

C*(X)[€0  +  €(X)  ^  €^(X)] 

=  C[€®  +  £(X)  -  £*(x)  ~  €^(X)],  X  e  Q  (14) 

where  Q  is  an  inclusion  or  a  subdomain  with  the  elastic  constant 
C*(x)  embedded  in  an  extended  body  D  with  the  isotropic 
elastic  constant  C.  €(x)  is  the  disturbance  strain,  €*(x)  the 
eigenstrain,  €''(x)  the  plastic  strain  and  uniform  strain  from 
which  O'®  =  Ce®  where  <r®  is  the  applied  far-field  stress.  The 
stress  O'  and  strain  €  are  vectors  expressed  as 

^  ~  0'22>  t^33»  O’ZS*  <3^13* 

€  =  {^119  €22,  €33,  €23,  €13,  €12 

where  notation  { •  •  •  represents  the  transposition.  Equation 
(14)  shows  that  the  eigenstrain  e*(x)  plays  the  same  role  as 
the  plastic  strain  e^(x)  does  in  the  homogeneous  materials. 

V^en  a  continuous  problem  is  discretized  for  the  purpose  of 
numerical  analysis,  the  domain  of  Q,  is  further  divided  into  N 
subdomains  ^,  (1  £  f  ^  N)  called  elements,  in  each  of  which 
boA  the  eigenstrain  c*  and  plastic  strain  are  assumed  to  be 
uniform.  The  stresses  in  element  ft,  caused  by  the  unit  eigens¬ 
train  €*  in  element  ft;  are  the  same  as  those  caused  by  the  unit 
plastic  strain  Thus  the  stress  o*-  in  ft,  induced  by  both  of 
the  eigenstrains  and  plastic  strains  in  all  the  elements  is 

0-1  =  1  Gije^  +  I  Gite?,  1  =s  j  s  (15) 

J  k 

where 


*  ^ 

^22*  €33,  €23*  €13»  <^12!^ 

==  €22,^33,  €23»  €b,  €a}J 


in  ft;  and  G,;  are  6  X  6  matrices  of  the  residual  stress  influence 
coefficients,  which  are  contracted  from  Gju„„,(ft, ,  ft;). 

For  the  generalized  plane  strain  problem,  W3  does  not  vary 
along  the  Xs-direction.  The  total  normal  strain  component  €33  is 
equal  to  zero.  The  strain  component  €  §3  is  here  taken  to  be  zero. 


Fig.  9  Conriparison  of  surface  plastic  strains  with  and  without  consider¬ 
ing  crystal  anisotropy 


Since  the  total  strain  el^  is  the  sum  of  e?3  and  £33,  hence  €33  s 
0.  The  stress-strain  relation  may  be  given  as 

erf  =  Cci  -  C(€f  +  ef)  (16) 

where 

"X  +  2/i  X  0  0  o' 

X  X  +  2^1  0  0  0 

7s^  X  X  0  0  0 

0  0  2fj.  0  0  ’ 

0  0  0  2fi  0 

-  0  0  0  0  2/u. 

'X  +  2p  X  X  000' 

X  \  +  2fj,  X  000 
r_  X  X  X  +  2j«000 
0  0  0  2m  0  0 

0  0  0  0  2m  0 

L  0  0  0  0  0  2m. 

and 

Ci  =  {fill.  *22.  *23.  *13.  *12  )f  With  €33  S  0 

From  Eqs.  (15)  and  ( 16),  expressing  the  disturbance  strains  ej 
in  ft,  in  terms  of  e*  and  we  obtain 

e,  =  C-[C(6.?'  +  €f)  +  I  G,ef  +  I  G*€j], 

j  k 

1  sisN  (17) 

where 

1  -V  -V  0  0  0  " 

J  -1/  1  -1/  0  0  0 

C"'=-  0  0  0  1+1/  0  0 

“000  0  l+vO 

LOGO  0  0  1  +  u. 

in  which  E  is  flie  Young’s  modulus. 

The  left-hand  side  of  Eq.  (14)  can  also  be  written  for  the 
generalized  plane  strain  problems  as 

<r,  =  <r*  +  C*€,  -  C*£f  (18) 

where  C*  is  a  6  X  6  anisotropic  elastic  stiffness  matrix  and  C* 
a  6  X  5  matrix  formed  from  matrix  C*  by  taking  the  third 
column  away.  And  erf  =  C*6®  in  which 

==  {^ii,  €221  ^23»  €13,  €12}^  with  £33  =  0 

Substituting  Eq.  (17)  into  Eq.  (18)  and  Eq.  (15),  Eq.  (14) 
becomes 

Te,?-  +  I  S,j€;  =  q,  +  (C*  -  T)€f  -  X 
J  k 

1  sisN  (19) 

where 

T  =  C*C‘'C 
S,  =  (C*C-'-I)G.y 
q,  =  a®  -  trT 
0-®  =  Ce® 

and  I  is  a  6  X  6  unit  matrix. 

It  can  be  seen  that  eigenstrains  «,*  depend  on  the  applied 
load  <r®  and  the  distribution  of  plastic  strains  e'  as  well.  Equa¬ 
tion  (19)  is  a  set  of  linear  algebraic  equations  with  unknowns 
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£,* .  Solving  these  equations  for  e,*  in  all  elements  fl, ,  we  can 
express  ef  as 

ef  =  P;  +  S  I  ^isN  (20) 

k 

where  p,  is  a  vector  in  terms  of  cr°  in  H/  and  R,*  is  the  matrix 
relating  the  plastic  strains  €?  in  Cly  to  the  eigenstrains  ef  in  fi,- . 
The  total  stress  in  Q/,  obtained  by  substituting  Eq.  (20)  into 
Eq.  (15),  is 

cr,-  =  cr®  +  <r  • 

=  +  I  G.yp,-  +  I  (Ga  +  I  G.yRy*)€?  (21 ) 

J  k  J 

In  the  above  equation,  the  term  of  2  GijPj  is  a  modification 
to  the  effect  of  the  applied  load  o’®  and  2  G,yR;*  the  one  to  the 
matrix  G,*  of  the  residual  stress  influence  coefficients  due  to 
the  inhomogeneity. 

Under  the  stress  component  alone,  Eq.  (21)  reduces  to 

Ci  =  5,0*22  +  X  GiJfeC* 
k 

where 

PioSz  =  +  £  Gi;P; 

J 

and 

G*  =  G/*  +  X  G,yRy* 

J 

For  slip  system  afi  in  element  Q; ,  we  have 
el  = 


forming  the  inhomogeneous  problem  into  a  homogeneous  one. 
The  summation  with  regard  to  slip  system  ^tj  includes  all  the 
slip  systems  in  fcth  element.  It  can  be  seen  that  Eq.  (22)  corre¬ 
sponds  to  Eq.  (9). 

5  Effect  of  Crystal  Elastic  Anisotropy  on  Fatigue 
Crack  Initiation  in  Polycrystals 

5.1  Elastic  Constants  of  Monocrystalline  and  Polycrys¬ 
talline  Crystals.  In  the  present  study,  nickel-base  superalloy 
polycrystals  of  FCC  structure  are  taken  for  the  numeric^  analy¬ 
sis.  Fatigue  crack  initiation  of  other  materials  can  be  analyzed 
by  the  same  approach. 

There  are  many  experimental  data  available  for  nickel-base 
superalloy  single  crystals  in  the  literature.  Since  the  alloy  com¬ 
positions  vary  with  different  materials  tested,  the  elastic  con¬ 
stants  obtained  may  slightly  differ  from  each  other.  In  the  fol¬ 
lowing  analysis,  the  elastic  constants  of  a  nickel-base  superalloy 
NijAl  single  crystal  given  by  Yang  ( 1985)  are  used.  The  elastic 
compliance  moduli  in  10~^/ GPa  with  respea  to  the  crystal  axes 
[001],  [010]  and  [1(X)]  in  the  same  notation  as  Yang’s  are 

Sii  =  Sun  =  1.01 

Si2  ~  Sn22  ~  “0.393 
S44  “  4Si2i2  “  0.848 

The  crystal  orientations  with  respect  to  the  loading  axes 
(fixed  coordinates)  are  represented  by  the  Euler  angles 
tp)  as  shown  in  Fig.  3.  The  global  elastic  compliance  in  tensor 
notation  is 


“COS  ^  cos  ^  sin  ^  -  sin  cos  ip, 
“COS  ^  sin  sun^  +  cos  ^  cos 
sin  ^  sin  ^ 


sin  0  cos  <f> 
sin  9  sin  <p 
cos  9 


and 

where  e^pi  and  Tapi  are  the  plastic  shear  strain  and  resolved 
shear  stress,  respectively,  and 

m^p  “  {fftii,  m22»  W233,  m23,  W13, 

“  {Wii,  ^22.  W33»  ^23*  ^13 »  ^12  }a^ 

in  which  m,y  and  riij  are  composed  of  direction  cosines  of  slip 
system  aP.  The  total  resolved  shear  stress  in  the  slip  system  is 

"^api  “  “T  affj  T  (g0!  *4"  T  apt  (2^) 

in  which 

Tagi  =  =  iafliOrIz 

=  S  Z  Cia0,  Q.j\  ^7), 

*  Cl? 

where  s^pi  is  the  modified  Schmid  factor  with  the  consideration 
of  anisotropic  elastic  constants  and  G(a/3,  fi,* ;  Qy)  expressed 
as 

are  the  modified  residual  stress  influence  coefficients  after  trans- 


^ijkm  “  ha^jP^k>JmySti0ky  C23  ) 

where 

The  inverse  of  Sijkm  is  the  elastic  stiffness  modulus.  Equation 
(23)  was  used  to  obtain  the  elastic  stiffliess  matrices  of  single 
crystals  referring  to  the  specimen  axes. 

A  range  of  test  data  of  Ae  Young’s  modulus  and  shear  modu¬ 
lus  of  NisAl  polycrystal  and  its  alloys  are  found  m  the  article 
by  Stoloff  ( 1989).  The  isotropic  elastic  constants  of  Poisson’s 
ratio,  V  =  0,3,  and  shear  modulus,  ^  =  65  GPa,  are  here  used 
for  our  analysis  of  NisAl  polycrystds. 

5.2  Physical  Model  and  Numerical  Results.  A  physical 
model  to  study  the  elastic  anisotropy  and  inhomogeneity  of  crystals 
on  the  fatigue  crack  initiation  in  a  polycrystal  is  shown  in  Fig.  4. 
An  aggregate  of  crystals  with  various  crystallographic  orientations 
is  embedded  near  the  free  surface  in  a  semi-infinite  elastic  medium. 
The  cross-sections  of  these  crystals  are  either  whole  or  part  of  a 
hexagon.  Since  the  crystals  at  some  distance  away  from  the  surface 
crystal,  in  which  a  fatigue  band  is  analyzed,  are  randomly  oriented, 
diese  suiTOunding  crystals  are  assumed  to  be  isotropic  and  homoge¬ 
neous.  The  fatigue  band  is  composed  of  three  thin  slices  P,  Q, 
and  R.  The  orientation  of  surface  crystal  #1  is  determined  by  the 
orientation  of  the  primary  slip  system,  which  is  45  deg  inclining 
to  X2-axis.  The  elastic  sti^ess  matrix  of  this  crystal  with  respect 
to  the  loading  axes  is 
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c*  = 


244.33 

118.19 

83.91 

-27.42 

-82.29  -13.71’ 

118.19 

244.33 

83.91 

82.29 

27.43  -13.71 

83.91 

83.91 

278.61 

-54.86 

54.86  27.43 

-13.71 

41.14 

-27.43 

153.56 

27.43  27.43 

GPa 

-41.14 

13.71 

27.43 

27.43 

153.56  -27.43 

.  —6.86 

-6.86 

13.71 

27.43 

-27.43  222.13. 

The  surface  crystd  with  the  fatigue  band  and  the  surroundins 
^t^s  are  shown  m  Fig.  4.  The  fatigue  band  is  divided  into 
paraUelogim  elements  and  the  remaining  parts  of  the  crystals 
are  divided  into  trapezoidal  elements.  The  paraleUogTam  ele- 
ments  are  assumed  to  have  linear  inelastic  strains  and  the  trape- 
roic^  elements  to  have  constant  inelastic  strains.  He.  5  shows 
the  layout  of  these  elements  in  the  aggregate. 

The  locd  stresses  in  die  surface  crystal  containing  the  fatigue 
band  are  affected  not  only  by  the  appUed  load  but  also  by  tte 
onenteuons  of  its  surrounding  grains.  Plastic  strains  can  develop 
both^ide  and  outside  the  fatigue  band  and  the  residual  stresses 
are  affected  by  *e  orientations  of  the  surrounding  crystals.  The 
globd  etetic  stiffness  matrix  of  ciystals  of  the  orientation  (9 
0,  ij/)  is  denoted  as  C*(9,  4>,  ^).  ’ 

The  surrounding  ^stals  can  have  different  sets  of  oriente- 
tions.  In  the  following,  one  set  is  calculated  to  illustrate  the 
ett^  of  muso^py  and  inhomogeneity  on  the  fatigue  crack 
imtianon.  In  this  set,  the  orientation  of  crystals  #2  and  #3  is 
chosen  to  be  (60,  30,  30  deg),  crystals  #4  and  #5  to  be  (0,  0 
O  d^)  and  crystal  #6  to  be  (30, 30, 60  deg).  The  global  elastic 
sufeess  matnees  of  crystals  wi*  these  orientations  are  listed 
as  follows 


C*(60,  30,  30  deg)  = 


the  sliding  starts  inside  the  band,  where  the  Schmid  factor  takes 
the  maximum  value.  The  case,  in  which  o-Jj  =  397MPa  r'  = 
±0.5 MPa  Md  T^=  200MPa,  was  analyzed.  The  surface  plastic 
resmved  shear  strain  vs.  the  number  of  cycles  is  shown  in  Hg 
7.  If  crystal  anisotropy  is  neglected,  the  fatigue  band  will  not 
TOvelop  in  the  surface  crystal  under  fte  applied  load.  Therefore 
the  toeshold  of  fatigue  crack  initiation  is  lowered  with  the 
considerauon  of  crystal  anisotropy.  Fig.  8  shows  the  plastic 
resolved  shear  strain  distributions  in  P  at  100  cycles  and  2000 
cycles,  respectively. 

The  cumulative  surface  plastic  shear  strains  versus  the 
number  of  loading  cycles  for  two  cases,  one  with  the  consid¬ 
eration  of  crystal  anisotropy  and  one  without,  were  calcu¬ 
lated  Md  plotted  in  Fig.  9.  It  is  seen  that  the  neglect  of  the 
crystal  anisotropy  in  a  polycrystal  may  cause  significant  er¬ 
ror  in  calculating  the  slip  distribution  and  in  estimating  the 
faugue  crack  initiation.  * 
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C*(0,  0, 0  deg)  = 
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C*(30,  30,  60  deg)  = 


247.36 
122.16 

76.91 
-9.89 
40.29 
L  3.64 
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-34.42 
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With  the  meth^  described,  the  modified  Schmid  factor 

Kofi  along  slice  P  is  shown  in 

rig.  6.  The  Schimd  factor  varies  from  0.42  at  the  free  surface 
to  the  maximum  value  of  0.59.  It  is  seen  that  the  effect  of  this 
^sotropy  on  the  response  to  the  applied  load  makes  the  Schmid 

hf  f  ^  free  surface  seen^  to 

constramt  on  the  deformation  near  the  surface 
f ‘“f  ^®”sidering  the  crystal  anisotropy, 

fte  Setod  factor  IS  umformly  0.5.  The  increase  of  the  Schrxtid 
factor  decrees  the  applied  load  required  to  activate  the  plastic 
slip  in  the  bMd  m  comparison  with  the  isotropic  and  homoge¬ 
neous  case.  The  region  in  which  the  total  resolved  shear  strL 
reaches  the  cnucal  shear  stress  will  slide  fim.  Therefore, 
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